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ABSTRACT 


Sound production and modification in an Argon discharge plasma 
was investigated theoretically and experimentally. Sound pro- 
duction by ions in the cathode region was shown to be feasible and 
consistent with earlier experimental results. A numerical analysis 
of the theoretical dispersion relation of Ingard and Schulz showed 
that plasma amplification effects will not compensate for acoustic 
losses. The theory of the mobility-limited thermionic diode 
operated as a microphone was expanded. Such microphones were used 
in an experimental investigation of sound production and modifi- 
cation by a discharge plasma. Possible plasma effects on sound 
propagation were observed. The characteristics of sound pro- 
duction by the plasma were found to be in reasonable agreement 
with theory. A departure from theory at low frequencies was inter- 
preted as arising from the effects of particle diffusion, an effect 


not previously observed in acoustics. 
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ee LN RODUGT TON 


Plasma acoustics deals with the interactions between charged 
particles and sound fields. The problems encountered are different 
from those of ordinary acoustics since charged particles experience 
bulk forces from external fields whereas neutral particles do not. 
When neutral particles are present, they are indirectly influenced 
by the external fields through their interactions with the charged 
particles. The complete mechanism of interaction is very complex, 
involving 10nization and recombination as well as elastic and non- 
feonizing inelastic collisions. 

This paper focuses on the problem by studying the effect of 
the interactions on the neutral particles. Two aspects are 
treated. One is the creation of sound by the current modulation 
eae discharge. The other is the modification of the pours when 
1t propagates through a discharge. The two aspects are closely 
related. 

The creation of scund in a gas requires a localized time-de- 
pendent pressure variation. Ina weakly ionized gas, the partial 
pressures of the ion and electron fluid components are negligibly 
small with respect to the partial pressure of the neutral gas. For 
this reason, the acoustic wave made in a discharge plasma is con- 


cerned only with the neutral component except when the effect of 


J : : ; 

"Sound! and "acoustic signal" or “acoustic wave! are used 
Synonomously and are defined to mean any fluctuation in the 
pressure of the neutral particles. 








ions and electrons on the neutral component is of interest. The 
pressure, expressed as p = KNT, Te be changed by either a change 
ip temperature, T, Or particie density, WN. 

In the positive column of a glow discharge, energy is con- 
mrnvally fed into the plasma from the electric field, which draws 
its energy from an external power supply. The energy is received 
by the electrons and ions in the plasma as kinetic energy. Since, 
in the aggregate, these particles are traveling at their terminal 
velocity, referred to as drift velocity, the increased kinetic 
energy of the individual electrons and ions is dissapated through 
collisions with the neutral particles. The averace rate of energy 


drawn from the field, per unit volume, can be written 
H =H +H, = eNU BE” + eN.U.B 
1 ee oe 


where om and Ht. are the electron and ion mobilities. Since the 
electron mobility is much larger than the fs mobility and the 
other terms are equal, the ion term ‘is negligible. Thus only the 
electrons need to be considered in the energy transfer process. 

In the equilibrium situation, the energy received by the 
electrons must be balanced by losses. The energy distribution of 
the electrons closely approximates a Maxwell Boltzman distribution. 
The difference is due to the effect of the electric field and the 
energy-dependent cross sections of the neutral particles with 
Whach the electrons interact. The electrons transfer energy to 
the neutrals through elastic and inelastic collisions. The energy 


transferred in inelastic collisions is usually transferred out of 
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the system of neutral particles through radiation. The energy 
received by the neutral particles in elastic collisions is thermal 
energy Which, at equilibrium is balanced by thermal losses to the 
container and neutral gas outside of the ionized region. This 
energy transfer between the hot electrons and cold neutrals is 
comparable to mixing two gases at different temperatures except 
that the hot gaseous component stays hot. The elastic collision 
losses are,normally much larger than the inelastic losses which 
can be neglected for this analysis. 

Thus the processes discussed suggest a mixture of two fluids 
of different temperature which retain their identity in the mixture. 
The hot fluid receives energy from an external source to balance 
meewconersy transferred within the mixture to the colder fiuid. the 
cold fluid is hotter than its surroundings, to which it loses 
energy. 

The characteristic times for the three different processes; 
energy gain by the electrons, transfer to the neutrals and dis- 
sapation to the surroundings, are different. The latter process 
is very slow. The characteristic time for energy transfer to the 
neutrals is of the order of the electron-neutral collision time, 
inomeally ~ 107° seconds. Thus, if the energy transfer to the 
electrons is caused to vary at acoustic frequencies, by varying 
N. or E, the time delay can be neglected and the energy transfer 
to the neutrals can be considered in phase with the energy transfer 
to the electrons. Since the time for energy transfer to the sur- 


roundings is much greater, the temperature change, and therefore 
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the pressure, of the neutrals will be essentially in phase with 
the energy transferred in. This means that a change in either E, 
mie weclectric field or No? the electron density, can generate a 
time dependent pressure component in the neutral gas which can act 
aS an acoustic source. 

Amplification of a wave occurs when the wave causes a change 
in the medium which leads to a change in the wa ve in phase with 
iselif, 1.€. "positive feedback" occurs. This can occur in a 
discharge plasma through the interaction of the pressure wave in 
the neutral gas with the electron SUE source term Hoe The 
neutral gas density changes in phase with the accustic pressure. 
Since both the mobility and electron density depend on the neutral 
density, the feedback can be introduced through them even with no 
change in the electric field. 

In the simplest analysis, the electron density varies in phase 


with the neutral density. If n, and n_ are the perturbations in 


the electron and neutral gas densities Ne and Ny? then 


This indicates that the energy transfer to the neutral gas varies 
in phase with the acoustic pressure variation. Since this would 
produce a pressure variation in phase with the original pressure 
variation, it constitutes an amplification mechanism. For net 

growth to occur, the acoustic losses due to thermal conductivity 


and viscosity must be less than the acoustic source gain. 








A more sophisticated analysis indicates that phase differences 
occur between the electron, ion and neutral density variations with 
resultant generation of a space charge electric field. The phases 
of all of these quantities influence the phase of the feedback and 
determine whether the plasma effects increase or decrease the net 


[mee TlAtLON Of the acoustic wave. 


A. CHRONOLOGY OF THE EXPERIMENT 

Interest in plasma acoustics here at the Naval Postgraduate 
School arose out of the study of moving striations. Earlier work 
performed here was that of Partlow tall Carreta and Moore EZ), 
Wilson ie Crandall [4] and Melville Bolg Melville's work 
represented a shift away from the study of moving striations. This 
was prompted by a theoretical paper by Ingard in 1966 [6], which 
predicted amplification of an acoustic wave ina discharge plasma. 
A more comprehensive theory by Ingard and Schulz followed in 
1967 L7]. , 

The experimental investigation described in this work began 
in 1968 as a natural continuation of the earlier work performed 
here. The earlier work here and elsewhere had demonstrated that 
sound could be generated by a plasma. This coupled with the 
prediction of effects on the sound as it passed through a discharge 
Medethe author to the consideration of a combined experiment. In 
the proposed experiment, a discharge would be used to create sound, 
and it would be observed before and after it propagated through the 
discharge. With the proper techniques, the experiment could provide 


PmLOrmation about both the creation of sound and its modification 


in the discharge. 


TS 








A numerical analysis of the Ingard-Schulz dispersion relation 
contained in [7] was performed as part of this investigation in 
which the variables were neutral density, neutral temperature, 
electron density, electron temperature and angular frequency of 
the acoustic signal. The attenuation was minimized with no con- 
straints placed on the values of the variables Seaeot for upper 
and lower limits. The results showed that the minimum attenuation 
occurred at high neutral densities, high electron densities and 
temperatures, low neutral temperatures and low frequencies. The 
maximum relative amplixication® was also investigated. The con- 
clusion drawn from a study of the behavior of these two quantties 


was that the experiment should test the theory in the medium 


observed by Yatsu1, et.al. [8,9] indicated the need to design the 
experiment so that the polarity of the discharge could be reversed. 
Also, it was felt that thermocouples were needed to determine the 


neutral gas temperature in and out of the discharge. 





ta 


"Relative amplification"(A_) is the change in attenuation 
due to the plasma amplification referred to the classical atten- 
uation. It is defined as ‘ 


where B is the attenuation coefficient, I Signifies ionized and u 
Signifies unionized. A negative value means that an amplification 
is Operative in the plasma case. 
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The experiment was designed to use two mobility-limited 
thermionic diodes as microphones to observe the signal at two 
points with the signal being produced by a remote source svegeeey 
Such an arrangement would permit a direct determination of the 
attenuation. The theory of the diode microphones was developed 
beyond the earlier elementary.form so that an extensive cali- 
bration was ae needed to use it as an analytical tool. Technical 
problems were solved so that the tube could be constructed with 
a smooth surface, free from side arms, bumps and bubbles which 
affect an acoustic signal, in spite of the numerous tube elements 
contained. A circuit was designed to control the discharge cur- 
rent to produce the sound in the desired form. 

molocercurrents@erc used to nroduce the seund with) the veource 
discharge. The current and sound waveforms were }ourler trans- 
formed to determine their spectra. By a process of comparison, 

Bees polation and comparison, the frequency dependence of the 
acoustic attenuation and the source transfer function were obtained. 
The measured values of the experimental parameters were used to 
determine the predicted values of attenuation from the Ingard- 
Schulz equations. These equations predicted that the plasma would 
not affect the attenuation under the conditions used. The experi- 
mental value for the attenuation co-efficient agreed with the pre- 


dicted value within experimental accuracy although the results 


indicated the possibility of a plasma effect. 








Lhest esults Lot Ya tsuumer nal [8,9] were then reevaluated. 
It was noted that while the relative change in the attenuation 
coefficient was small, since the attenuation was very large, the 
amplification effect was large. In view of these results, it 
was recognized that the quantity to optimize is the difference 
in attenuation co-efficients between the ionized and un-ionized 
cases. 

Earlier a theoretical study of the frequency characteristics 
of a discharge sound source was undertaken by Schulz Lio]. He 
concluded that the energy transfer from the electrons was the 
dominant source mechanism. The predicted frequency dependence 
is sin (vta/c)/(vta/c). When the source transfer function was 
experimentally detetmined, it was found to be consistent with 
theory at high frequencies but not at low frequencies. The low 
frequency effect was attributed to diffusion effects. 

in the following two sections of this haste earlier experi- 
ments relating to the creation and modification of sound are 
discussed. Following, Chapter III discusses the theoretical treat- 
ment of sound creation and modification. This order of presentation 
was chosen because it was felt to offer a better insight into the 


aspects which the theory should treat. 


B. PREVIOUS EXPERIMENTAL WORK ON ACOUSTIC INTERACTIONS 
Before considering specific papers of direct interest to this 
experiment, a quick general review of previous experimental investi- 


gations of the acoustic interaction is in order. Interest in the 
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interaction between charged particles and neutral acoustic waves 
is recent. A relationship between moving striations and acoustic 
waves was suggested by Watanabe and Oleson C11] abo mel SI ova re 
Wojaczeck iien reported observations of a sound interaction in 
1960 during experiments in Xenon. He noted on effect in his main 
discharge which occurred when an auxiliary discharge was pulsed, 
and he identified the effect as due to sound waves generated by 
the pulsed discharge. | 
Since that time, many workers have reported observations of 
sound waves generated by discharges. The source and type of 
disturbance in the discharges, which were observed to produce 


sound, have been sufficiently varied to conclude that it is a 


1) 


general phenomenon. Sound has been observed Li-3] whose source 
was interpreted to be the current fluctuations associated with 
Poeveng Striations. Sond production by artificially introduced 
current fluctuations has been reported Pe 15 for a variety of 
conditions. Modulation of the electron temperature in a discharge 
by an R.F. field was observed by Fitaire and Mantei [16] to produce 
sound. Multiple pulse emission from an HCN laser, reported by 
iheimer, et.el., eile was interpreted by Born [18] as due to sound 
production by the discharge. 

The variety of techniques used to study the sound produced by 
pulsed discharges is almost as interesting as the observations 
reported Pies o=2 > |. The pulses ranged in duration from one to 


: -2 
fifty microseconds [23,18], in pressures ranging from 5 x 10 torrl19] 
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Oe cO, VOLE eaesele In some experiments, the sound was detected by 
observing its effect on some property of the discharge which 
produced the sound. An example of this is the work reported by 
Berlande, et.al. ee in which they observed perturbations in 
the light of the afterglow of the plasma which produced the 
sound. Nygaard [21] utilized. the light perturbation effect by 
uSing one Aeracae to produce the sound while he observed the 
afterglow of a second discharge. Goldstein eed, Les) weed 2 
mobility-limited thermionic diode as a microphone. These and 
other techniques were used in some cases because of the problems 
associated with conventional acoustic detection techniques under 


the conditions involved. In other cases, they wece used because 


ct 


Mec Coren ts detect On Mia e anmiee so Sctod mecult 2m 


experiment. 

Only a few papers nave been published concerning the effect 
ofa Pp sma OnesOund pr Opagating through 1: [5,8,9,261, as 
opposed to the effect on the discharge. Because of their relevance 
and sparsity, these papers will be discussed in detail in the next 


section. 


See or eCIFIC EXPERIMENTS OF INTEREST 

The two aspects of this experiment are the creation of sound 
and its modification in a gas discharge. Earlier experiments have 
examined both of these and as the experimental results have been 
reported, theories have been developed to explain the experimental 
results. The theoretical aspects of the experiment will be dis- 


cussed in the next chapter. 
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1. Sound Production 

A review of the existing experimental evidence concerning 
sound production can point out features which must be explained 
by any theoretical model of the process. 

Wilson L3] determined from phase measurements of the sound 
produced by moving striations, that the only sound source position 
consistent through the experiment was at or near the cathode. By 
comparing the current and sound waveforms, he concluded that the 
Setnd was produced by fluctuations in heating processes occurring 
near the cathode due to current fluctuations. He suggested that 
the mechanism must involve an interaction between the fluctuations 
in ion density ANE temperature and the neutrals. 

Nygaara and Meitz LZij ace more Specific Aas tO the Location 
Srethe source. They state that the acoustic disturbance from their 
H.V. pulse discharge originates in the negative glow region. They 
postulate a high pressure driver region near the cathode maintained 
by electron heating and momentum transfer from the ions. 

Hayess [20] also interpreted the sound production as arising 
from localized heating of the neutral gas. However his experiment 
brings in a different aspect. He used a high voltage pulse applied 
to an electrode near a discharge already running, so that the other 
electrode of the circuit had a much greater charge carrier con- 
centration than the exciting electrode. He found that the sound 
source was always at the exciting electrode, regardless of the 


polarity of the voltage pulse. By varying the rise time of the 


Ly 








voltage pulse, Hayess established that it had particular importance 
in the sound excitation process. Pulses with rise times of 

5) aoe 10” V/sec produced sound, those with rise times of 1 eae 10°v/sec 
did not. In addition, Sine waves with an average rise time of 

ex 10/ V/sec did not lead to sound production. Hayess also 
established that the presence of charge carriers at the exciting 
electrode reduced the sound produced. From consideration of these 
factors, Hayess concluded that energy gained in the high field by 
the charge carriers is transferred to the neutrals. He further 
concluded that the ions were responsible for the inertial phenomena, 
i.e. dependence on rise time. 

In contrast to the above experiments, the observation of 
sound nroduction by modulation of the electron temperature ina 
discharge, reported by Fitaire and Mantel [16], may be deemed as 
experimental confirmation of the roll of the electron temperature 
in the acoustic source term of Ingard and Schulz's equations AIF 
In interpreting these experimental results, it must be borne in 
mind that there is a distinction between the energy of electrons 
in a low field region, such as the positive column and a high 
field region such as the cathode fall. In the latter, the motion 
is highly directional and can not be considered thermal. 

Upon considering the above experimental results, the 
following conclusions inay be drawn. There are two mechanisms 
through which sound is produced in a discharge. One mechanism 
occurs near the cathode and probably involves ions in a high field, 


the cathode fall. Nygard's contention that the source is the 
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negative glow region is in contradiction with the results of 
Hayess concerning the high field. The other mechanism occurs in 
the positive column and involves the electron temperature. The 
former mechanism appears to dominate the production process. 
Further comment shall be reserved for a discussion of sound pro- 
Guctizon in Chapter II. 

ee SOUuld HModitlcatlon 

The experimental evidence conte no the effect of the 
plasma on sound is not so conclusive. Wojaczek [12] observed the 
Merturbation, produced by the sound, in the light emitted from a 
Xenon gas discharge. On the basis of his calculations, he con- 
cluded that the attenuation was much greater than expected. He 
attributed tnis to pussibie interaction with ¢ 

Later Hayess [19] pointed out that the attenuation 
Wojaczek predicted was in error because the dominant loss 
mechanism, wall effects, was neglected. avecemeeconpated the 
predicted values. With the wall losses taken into consideration, 
the experimental values agreed with the predicted values within 
the experimental uncertainty. Wojaczek's experiment was. conducted 
fameressure from .0O7 torr to 1.5 torr with currents from 1.5 to 
O amperes. 

Working in Argon, Yatsul1 et.al. [8,9] studied the effects 
of a discharge on ultrasonic sound. They used barium titanate 
transducers with resonant frequencies of 50 kHz, 100 kHz and 
200 kHz for transmitters and receivers. The signal received with 


no discharge, as a function of transducer input voltage, was the 
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reference. By dividing this into the results obtained with the 
discharge on,a "spatial growth rate' was obtained. A definite 
change in the growth rate was observed. It was found to be de- 
pendent on the current and pressure of the discharge and on the 
polarity with respect to the propagation direction. The polarity 
effect was so strong ee under conditions such that while propagation 
in the direction of electron drift motion enhanced the signal by 
50%; no signal was received when propagated contra to the electron 
drift. The variation of the growth rate with pressure and current 
was not monotonic, as the signal exhibited a maximum enhancement 
at intermediate values of the range investigated. The conclusion 
drawn by the authors was that strong damping occurred for contra- 
pronagating stanals under all conditions. For parallel propvaaation. 
there 1S a critical pressure pe above which attenuation occurs and 
below which amplification occurs. They attributed the directional 
effect to coupling between the sound and A ne motion of the 
electrons, a coupling which is strongest when the sound velocity 
1s nearly equal to the electron drift velocity. 

Their experiment provides experimental support to the con- 
cept of energy transfer from the electric field to the neutrals 
via the electrons. It should be pointed out that the acoustic 
frequency is of the order of the collision frequency at the low end 
of their pressure range, hore to .5 torr. Indeed it was at this 
end that the maximum effect was observed. 

Macon [26] used a Helium discharge to investigate the 


effect of the discharge on sound. He employed resonant cavity 
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techniques with frequencies below 2000 Hz. The sound source was 
calibrated to provide constant amplitude deflections over a range 
of frequencies. A rigidly mounted microphone received the signal 
at the other end of the tube. By comparing the ratio of pressures 
at adjacent pressure maxima and minima, and assuming the net loss 
was dominated by boundary effects, the attenuation constant was 
determined. The value for the "discharge off'' condition was com- 
pared with the "discharge on" condition. Measurements for the 
discharge off condition were taken one-second after turning the 
discharge off in order to have comparable temperature. The tem- 
peratures, as determined from sound velocity measurements, dropped 
from 450°K with the discharge on to 360°K with the discharge off, 
with one-second. The results indicated a decrease in atteuuation 
rate of 30%; a value too large to be accounted for by the pressure 
and temperature differences. Although the experimental results for 
attenuation and amplification rates were both in poor agreement 
with theoretical predictions, Macon concluded that amplification by 
the electrons was definitely occurring. 

These experiments indicate that there are amplification and 
attenuation effects. Wojaczek's experiment does not hint of any 
Srirectional effect. Macon did not look for any. The directional 
effect observed by Yatsul is explainable by models for other fields, 
namely electronics and solid state. The Ingard-Schulz model to be 
discussed in the next section does not take directional effects 


mito consideration. 
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Physical theories go hand in hand with experiments. Theories 
are developed to account for experimental results. Experiments 
are conducted to evaluate theories by testing their predictions. 
In the previous chapter, the experimental evidence regarding pro- 
Suetion and modification of sound by a discharge plasma was re- 


viewed. [n this chapter, the theories will be examined. 


A. SOUND PRODUCTION BY DISCHARGE PLASMAS 

After several investigators reported observation of sound 
produced by electrically modulated discharges, Ingard [6] per formed 
a theoretical analysis of the problem. He derived the wave equation 
Hor the orossure in the neutral gas in which the interactions with 
the electrons and ions were expressed in terms of sources in the 
wave equation. Neglecting the fluctuations in neutral particle 
Pefesty resulting from the imbalance of ionizaticn and recom- 


bination rates, Ingard obtained the following wave equation: 


: . 
10°» 2. _ y-1 OH 7 
Boye Ue > ae ean (2.1) 


In this equation, p is the acoustic pressure, defined as the time- 
dependent pressure variation about the ambient pressure; c, the 
Speed of sound; and Y is the ratio of specific heats. "H' is the 
energy source term, representing the rate of energy transfer to the 
neutral gas from the charged particles per unit volume. The rate 
of transfer of momentum to the neutral gas per unit volume is re- 


—_ 


presented by the momentum source term F. 
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Ingard considered the effect of momentum transfer represented 
by V-F as negligible for two reasons. First the electron and ion 
contributions tend to cancel each other. Second, only the spatial 


—> 
variation of F leads to sound, and for most cases, this is very 


Small. For the energy source term, he wrote 
Am mv | 
H = —— | —-——/ nN (oO) (2.2) 
mn Z eon 


iopere we is the average electron velocity and on a TEES average 
eellision cross section for elastic scattering. Ingard considered 
the 10n contribution negligible due to its much lower mean energy. 
Ingard pointed out that since the electrons were in an equilibrium 
state, the energy lost to the neutrals, represented by H, was the 


Same as the energy gained from the field. Therefore 
H = N cE (4 JE) (2eno,) 

where Uo is the electron mobility and E the pileciaie sete Jual 

Schulz [10] used equation (2.1)*’as a basis to develop an 
expression for the sound field produced by current-modulating a 
discharge operating between two planar electrodes. He considered 
the effects of the two source terms separately in calculating the 
pressure response. For both, a one dimensional problem was 
assumed and it was assumed that the time dependence was given by 
exp(-1Wt). Under these assumptions, the wave equation inside the 
Source region becomes 


ae ) ‘ Gd | 
3 (2) po =a 3 (Y-1)H. (2 24) 
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If H is assumed to be spatially constant within the region, 


then the general form of the solution inside of the region is 


~ ik i ae 
Pe cman at oS + i poo Ble (25.5) 
in UW 


With the additional assumption of symmetry about the origin, 
chosen to lie between the two electrodes which are separated by 
a distance a, Schulz's expression for the pressure field outside 


mae Source region for the energy source is obtained. 


-1 ; CLA ; ob 
P = =). H sin te exp(ilz|%). (2cy) 
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Senuiz assumed tha= the momentum transrer acoustic Source 
meese from an unbalanced fluctuation between the forces exerted 
by the ions and the electrons on the neutral gas. He stated that 


this force per unit volume was given by 


ee 2 Ww E 
FO, o 
tye W  -1W 
en 


: (2.7) 


With the assumptions W, ?7W and Won Yin? this simplifies to 


po iN. gE (W/W, _) ; (ZnS) 


~~ 


imthne electric field 13 1S assumed to be constant between the 


two electrodes and zero outside, then V°F can be expressed as 


V-F 


F [6 (245) : 8 (2-5) 


where 6(z) 1s the Dirac delta function. This leads to a pressure 


freid of : 


mca UA Z 
Pe =o wl to B sin (2) exp (t 1) 4 , (21) 
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Schulz concluded that the momentum source term would be negligible. 


He compared the two sound fields and found that the ratio was of 
Corea \ | 


3 — in favor of the energy transfer process. 
n / 





the order of 


Schultz's results predict that the energy transfer mechanism 
dominates. This result is in agreement with Ingard's earlier [6] 
Senclusion. Schulz's ec Sian for the sound field predicts 
that the frequency dependence is like ane where x = Vv — : 
Heretofore this frequency dependence has not been experimentally 
investigated. However there is experimental evidence which refutes 
Schulz's results and contradicts Ingard's assumptions concerning 
the energy source term. 

; the experimental evidence suggests that 
the principal source term lies near the cathode and involves the 
field. Schulz's analysis implied that the positive column was 
the source region. The nS frequency dependence also shows 
that the sound produced by high frequencies will be small. The 
effective frequencies of the sound pulses mentioned earlier have 
been in the range of 20 kHz to 1 MHz, all of which are high 
acoustic frequencies. Perhaps the comparison with Bee dis- 
charges is not a valid one because of the difference between the 
experiments and the conditions assumed for the derivation. The 
derivation presumes small perturbations yet the current changes 
in a pulsed breakdown discharge can not be small perturbations. 


Since sound production appears to occur in the vicinity of the 


cathode and appears to be associated with high fields, a possible 
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source term compatible with experimental results would be a field- 
dependent energy source term. A field-dependent source 1S one in 
which the energy for the acoustic wave comes from the field. The 
dependence on this term ee through the energy equation in 
contrast to Schulz's treatment, in which the field entered through 
the momentum equation. 

When the processes occurring in the region of the cathode are 
considered, a mechanism for the production of sound suggests 
itself to the author. The contribution of ions to the energy 
transfer is usually neglected. Due to their higher mobility the 


electrons transfer energy from the fieid to the neutral particles 


much more readily. It should be recognized that due to their 
Meemen macses, the ions and neutrals wall be at essentially Spe 
same temperature. The unity mass ratio means that during a 


Collision, an ion will transfer a large fraction of its energy. 

In the cathode region, this factor becomes important. There the 
10ns are the majority carriers of the discharge current. Their 
density is higher and they accelerate toward the cathode under the 
high field of the cathode fall. The increased energy will be 
teansferred to the neutrals through collisions. There will also 
be a net momentum transfer since the ions acquire momentum from 
jhe field and transfer it also. 

A rough approximation for the strength of such a source term 
can eaSily be developed. The average energy gained from the field 
between collisions is eEl. | where E is the electric field strength 
and L.. 1s the mean free path of the ions. Since approximately 


half of the ions energy would be transferred ina collision, the 
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energy transfer rate per unit volume is 


H. =~ + ch, W. Na. | (2.10) 
a Z pe ges 9 as 


This approximation is low because it uses the average energy that 
mould be transferred if the ions were stopped after a collision. 
It does provide an order of magnitude estimate. If oe is 


approximated by 0/0, then 


Dep mee Oe (2 las) 


The ratio of this term, where E Ba the field in the cathode 
fall to the electron source term, E = Ey? the field in the positive 
column is 


H, eB, viN, 
oe (2212) 


al 
H. N oles Meo eS) 
Now ee. is the cuxrent density in the positive column, and 


e EN. 1S approximately the current density in the cathode region. 


Since the two current densities should be the same, 


i 1 Be 
H 2E (ky 
€ p 


The approximations used in arriving at this result are equivalent 
Gomsaying that in the positive column, all of the power delivered 
goes into heating the neutral gas and in the cathode region, part 
goes into heating the neutral gas. The ions accelerate under the 


menron of the field in the cathode fall and gain kinetic energy. 
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The difference between the energy taken from the field and the 
kinetic energy gained represents the energy transferred to the 
neutrals through collisions. Since the current density is con- 
stant but the field varies throughout the length of the discharge, 
the energy transfer rate 1s a function of the field and the 


fraction of energy lost in collisions. 


B. SOUND MODIFICATION BY THE DISCHARGE 

The theoretical problem of acoustic waves in a weakly ionized 
gas has been studied by several authors other than Ingard eee le 
Without exception, the work of these authors needs modification 
before it can be used to predict the results of an experimental 
study ina laboratory. This is because they neglect boundary loss 
effects; yet at the low pressures of most discharge experiments 
with typical tube diameters of less than, or of the order of, 5 cn, 
boundary losses are dominant. The work of Ingard and Schulz bea 
is the most applicable to an experiment. In addition, it is a 
general treatment and is reasonably susceptible to numerical 
evaluation. 

The mechanism of amplification is embodied in the concept of 
energy transfer from the electrons to the neutrals through col- 
lisions. In the absence of an acoustic wave, the energy transferred 
is spatially uniform and merely serves to increase the temperature 
of the neutral gas. In the presence of an acoustic wave, the energy 
transfer becomes time dependent and spatially non-uniform. This 
produces an acoustic wave which can be in phase with the initial 


wave, leading to amplification. 
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Ingard's initial analysis of the problem L6] followed the 
experimental observations of sountwonoatene by electrically 
modulated discharges Pee ee ani He treated the problems as 
belonging to two categories, sound produced by current fluctuations, 
or an acoustic wave instability. For the first category, he derived 
equation (2.1). To treat the second category, he placed a term 
representing the electron density perennation into the energy 
source term of equation (2.1). His results indicated that the 
acoustic Signal from a current fluctuation and also the growth 
of an acoustic fluctuation would be small at room temperatures but 
enhanced at cryogenic temperatures. 

Although pupae ned in May 1966, the development of reference 6 
was done in early 1965. in the rail of 1966, macon 126} compieied 
his investigation. A thesis student of Ingard's, he attempted to 
conduct his experiment at cryogenic temperatures but had to use 
room temperatures because of noise due to moving striations. He 
concluded that amplification was falas place at temperatures 
greatly above room temperature (450°K). His values of growth and 
attenuation rates were not in agreement with theocy but he attri- 
buted this to his experimental arrangement. The more comprehensive 
treatment of the problem by Ingard and Schulz [7] was published in 
1967. 

The Ingard-Schulz analysis is based on a three fluid model of 
the ionized gas with both energy and momentum transferred between 
mae £Lluid components. The effects of Seeoene and heat conduction 


within each component and with the walls are also taken into 


Qo 
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consideration. The conditions assumed in making approximations 
may be considered typical for a glow discharge: gas pressures of 


ae O 2 - 
1-10 torr, electron densities of TO - iO) cm ae and acoustic 


frequencies 105 Ws5s 10°. 
The equations forming the basis for the analysis are introduced 
as linearized space-time Fourier transforms of the continuity 
equation (2.14), the momentum transfer equation (2.15), the energy 
balance equation (2.16) and the equation of state (2.17) for each 


fluid component (these equations are presented here for complete- 


ness and all items will not be defined at this time). 


- iun. + ikN.v, = 0 (2.14) 
j jj 
ee = ei ea mea, ws. (cmt 
~ —— >. = (een i -v. iene) << wyy. (2.15 
2 ea ee SS ao job 5 ‘es (a 
(j,k,£ = n,i,e and permutations) 


ee .5°5 = Ce Nee Baca, Sa) a 
- Fa ;(k We 0: Genet) : Wet 7) 
=) = - ik(KT/m.)v. - F5(k,) oa iOn 
=} 
Wee [CP Pn*P iP i)%o/Te i oh OG es miP (4%) | 
=, OP oP Bi 
a) aia —— 9.=c. KO. : 
= (3) = mans* (3 o 3 5 Y; a uae (22) 


(altee) 
The frequency W is understood to be real and positive with the 


waveform taken as exp (ikx-iwt). This means that a sign difference 








between the real and imaginary components 


constant k corresponds to amplification. 


of the wave propagation 


Amplification occurs when 


the amplification mechanism dominates and net growth of the wave 


Occurs. 


In the momentum transfer equations (2.15), the term $7505 0)y 


represents the viscous stress.in the j-th fluid component. A 


familar term f..(k,W)@. accounts for heat conduction in the j-th 
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component, where ie is the temperature perturbation. Poisson's 


equation is used to relate the electric field perturbation to the 


density perturbation. 


Ingard explicity assumed that inelastic collisions and ion- 


neutral energy transfer could be ignored. 


Other assumptions are 


iiomtert tit tice ANAL¥SIS Because 
mbectric field 1s not present in 


ambipolar diffusion toc the walls 


ef lected terms; the d.c. 
the energy balance equation, and 


is neglected. 


By algebraic methods, the equations were reduced to a set of 


three homogeneous equations for the component velocity perturbations, 


Pirteten here in matrix form. 


A B Cc Vv 


(2.19) 


Approximate expressions for these co-efficients were given by 


Meeard and Schulz, in the main body of their paper, as being 


Sy, 








consistent with the acoustic mode of motion in a gas with a Low 

















degree of 1onization. These equations Ader 
ors : a) es ne - 
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Weing the co-efficients given in equation (2.20), Ingard and 
Schulz made additional approximations and expressed the dispersion 


mPelation in the form; 
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It was stressed that this expression was valid only if the 
second and third terms in the brackets were small with respect 


poeunity. 








In the above expressions, the adiabatic speed of sound for 

the j-th component, os TSsae fi medsaas: 
2 "2 
Ce (02 IE) ; 

where: Y; = C5 Geer the ratio of specific heats; K, the Boltzman 
Sonstant -; We and We are the temperatures and mass of the j-th 
component. The term Ores is a constant of proportionality 
which is related to the effective collision frequency of a 
particle of the j-th type moving through the fluid component 4% 


with Ny particles per unit volume. 


The Kirchoff wave attenuation caused by viscosity and heat 





conduction, which was introduced through the terms fon and ie 
MEME M@cutral Gas 15 taken into Account by usc cf ti characteristic 
frequency Wn? defined by 
WD i 1 Weak 
See = 4 D220 
Ww U) T Ww is V2 Joy) ( ) 


The chatacteristic frequency for ions, Was is defined ina 
Similar manner. The ion and electron plasma frequencies are 
denoted by WW and Oe The term 1/WT expresses the effect of the 
heating of the neutrals by the electrons. 

im the analysis of their results, Ingard and Schulz point out 
that wave amplification will occur if the term in the imaginary 
Bart of eo Eq. (2.21), involving T dominates. They give a simple 
argument and develop a criterion for amplification involving the 


electron-neutral collision frequency Wan’ 
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They also point out that contrary to the result obtained by 
Ingard [6] earlier, this more general treatment indicates that a 
decrease in neutral gas temperature does not always ee: the 
amplification. They show that the optimum temperature can be close 
to room temperature, based on the approximate expression Eq. (2.20). 

Solving equation (2.19) for the ion and electron perturbation 
velocities in eae of the neutral velocity and using approximations, 
it is shown that the assumption often used, that One Ue and vi, ore 
all the same in the acoustic mode, is not generally valid. [t is 
true only for frequencies below some characteristic frequency. At 
higher frequencies, the ion and electron velocities decrease in 
amplitude and are brought out of phase with the neutrals. 

Macens observation of an amnlification effect at increased 
temperatures does not contradict the Ingard-Schulz development; 
@ienobservations of Yatsui et.al. do. Since Ingard-Schulz did 
not Boneider directional effects, and did not set up the equations 
so as to include any drift motion, their results could not have 
predicted Yatsui's results. The question apace whether the 
directional effects are present at lower frequencies and higher 
pressures. In the Yatsul experiment, the frequencies were high 
and the pressures low, resulting in conditions where the wave- 


lengths were on the order of the mean free paths. 
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IIIT. NUMERICAL EVALUATION OF INGARD=-SCHULZ EQUATIONS 


The deSpersion relations developed by Ingard and Schulz predict 
a change in attenuation, possibly even amplification for an acoustic 
wave in a weakly ionized gas. To improve the chances of observing 
what might be a small effect, it was felt that the optimum con- 
ditions must be determined and the experiment conducted under 
Similar conditions. "Conditions" here refers to the values of 
ccontrolable discharge parameters. 

The purpose of the experiment was to measure a change in the 
attenuation constant. When the terms in inGgardeaneemcelmbt zou 
expressions were analyzed, the obvious variables were the micro- 
scopic narameters: nentral narticle density ae electron density Na? 
neutral temperature Ty and electron temperature a and the angular 
frequency of the acoustic wave W. The problem of determining the 
optimum conditions was thus the mathematical problem of minimizing 
a general function of five variables with the variables constrained 
to a certain region, the region of experimental feasibility. This 
1s best done with a computer. 

Since as presented, the Ingard-Schulz equations were not 
Suited to such an analysis, the process was one of analysis and 
Synthesis. Analysis was needed to express the terms as relations 
involving only the above variables and physical constants of the 
gas such as cross section or molecular diameter. These terms were 
then evaluated numerically and Sohpined with each other to compute 
the predicted value of the attenuation constant for the input 


values of the variables. 
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Using a computer, the equations were evaluated for a wide 
range of values of the input variables. This provided a more 
thorough evaluation of the Ingard-Schulz equations than they had 
made. The results showed that their approximate expression, 
equation (2.21) predicted net amplification over a large region of 
variable space although the more complete Bnecenions never did, 
even when Ingard's criterion for the validity of the approximation 
was satisfied. 

The process of analysis is described in Section IIIA following. 
The procedure used in evaluating the terms and determing the 
Sereinum values for the variables 1S described in section IIIB. The 
results of the evaluation are presented and discussed in section 
ie. 

It should be pointed out that this section discusses the 
evaluation of the Ingatd-Schulz equations before the experiment 
was conducted. Afterwards, the program was used to compare the 
predicted values of the attenuation constant with the observed 


values by substituting known or estimated values for the variables. 


A. ANALYSIS 


1. Obtaining the Attenuation Factor 


The system of homogenous equations (2.19) has a non-trivial 
Selution Only if the co-efficient matrix is Singular, i.e. its 


Geterminant iS zero esas We can therefore write: 


|co-efficients|| = A(E1-HF)+B(GF-DI)+C(DH-GE) = 0 (ely 
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We define the complex quantity n = ck/W. It can be considered 
opr : : : 2 
an acoustic index of refraction. Since in Eq.(2.20) n appears 


only in the co-efficient A, we may solve for it by defining Ay as: 





Z a ae 
= = +7 foe 
Ay Atn i: 7g 5 
nn 
and substituting this into equation (3.1). After a rearrangement 


coe 
Of terms, n is found to be: 


nf Shek: CUDH=EG) - BUGE=Dr) (2) Dy 


I [1S sist EL-HF 


Now ao Ganpe written in the form 
Zs ; 
n = \+iv (ue xeall) 


and the desired result is the wave propagations constant k in the 


from: 
k = &+iB, (2,8 real) C3230) 


so that the problem reduces to solving for the square root of a 
complex number. In particular, since @ is defined to be positive, 
the root of interest is the one for which 8B has the same Sign as bb. 
This means that if ae is written in polar form: 


2 2; 19 
Nees all te le 


3 
mae solution is the root at 0/2, i.e. the principal root. The 


attenuation factor, 8B, the ultimately desired result, is given by 


B = = im(n,) (3.4) 


where ny ise the prineipaimuoot of ae. 


39 








2  ACOUSi eros se> 

An acoustic wave is damped through two physically different 
mechanisms, viscosity and thermal conductivity. The damping can 
also be distinguished by the location where the losses occur, 1.e. 
at the boundaries or in the bulk of the gas. The latter charac- 
terization has the advantage that it divides the losses into two 
terms of different frequency dependence. nee the losseagarc indian 
the two effects are not simply additive but instead are coupled 
together in the differential equations for an acoustic wave. As 
mers would only be the case for viscous fluids or for tubes of 
small diameter, neither of which are treated here, we shall con- 
sider the Becen be additive. Since the mean free path and 
therefore the boundary layer is inversely proportional to density, 
at sufficiently low pressures the assumption of additive effects 
may have to be re-examined, 

A plane wave traveling ina duct attenuates expontentially 
with distance with attenuation aoe hen: "at. With small losses, 
a is the sum of the buik attenuation co-efficient a_ and the wall 


B 


attenuation co-efficient ay: For a Civcular pipes these coq 


efficients are given by Morse and Ingard [35,p989] AS. 
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where H, the thermal conductivity, is given by: 


2D 
Soe aa ee 
and 
WS ee is the co-efficient of viscosity. 


VY 


The mean free path & is 


oe Ces 


2 nen 
In Riese eee Oe ioe thee mase sGens 2ey se ethemadi a (tlem peed 
of sound and D the atonic diameter. When these expressions are 
used for a three component fluid, c will be written C3 however 
there is no difference between the two, Co. =, 

Since the energy of a wave is given by the square of its 
amplitude, the energy of a wave attenuates with co-efficient 2a. 
Wai a velocity of a, the characteristic time for energy dis- 
Sipation by the neutral component, Lier is related to the atten- 


uation co-efficient a by: 


lees 
7 Ta a e (367) 
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Therefore, for the small damping approximation of a plane wave 

in a circular pipe, the characteristic time for energy dissipation 
: : 3. 

mengiven by”: 


me (AIA nm 2 ae 2 
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Sthe expression does not exactly agree with the expression in 
Ref. 7 but it is consistent with the derivation there and agrees 
with the same expression in Morse and Ingard [35]. The difference 
may be due to an error in typesetting. 
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Substituting for K and 7, using Y = 5/3 (a good approximation for 


a noble gas), the above expression simplifies to: 


ae 2h Si yaie [2-4 
F = 2-4 w eer || 2-4 Whe. : (3.9) 


This term enters the dispersion relation as W/W defined by 





Pemation (2.22), 








cee  ZCLrO Ionization Limit 
In order to be consistent, the attenuation predicted by 
migard’s equations, for the limiting case of zero ionization, 
should be comparable with previous results for acoustic atten- 
uation. For the standard of comparison we snall use the 
attenuation co-efficient "a" which was used as an input in these 


equations. 


29) B 
U us 5 L 
a a) [4 2 3 * be (2.4 “eS ; (3 0)) 
nN) 


Bole clomlonization, equattonas(2,21) ance s92 8 rednce 
to the same expression for nee 


ne 2 @ eee. (3.11) 


Now, both equations (2.21) and (3.2) were based on 
approximations. When the expressions for the co-efficients, 
Eq (2.20), were presented, Ingard and Schulz had already made 


Several approximations. Going back to the full expressions for 








the co-efficients, which are contained in the appendix of Ingard 
and Schulz's paper, an examination shows that when the ionization 


goes to zero, only one independent co-efficient Ay remains, given 


Dy : 5 
fa oe os 1- y-1 : ee ee te 
fe) Ww ye bse at fof Ww - 


Even though there is only one co-efficient, it must equal 


zero for the equations to have a non-trivial solution for V_. 


k2c2 " 
Setting ae equal to zero, and solving for n* a aa , we obtain 
Ww 
ae 
5 iti fy 7 
n Same is ; (orad2 ) 
Ve Vie eo) 
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where Beer and ie are the viscous and thermal conduction com- 
ponents of the acoustic losses given earlier in section iTA. 

They are hidden in the expression for the attenuation co-efficient. 
The attenuation time ie can be expressed in terms of these com- 
pOnents as 

-1 
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Using the same expressions and approximations as before 


they can be written as: 


Pee ee 
ae 32 St E\ [2-4 wee 
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Since there were fewer approximations involved in 
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I 


Eq. (3.12), the value of B obtained is a strongex comparison test 


emer the use Of equation (3.11). Using both equations (3.11) and 
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(3.12), 8 was calculated for several sds fer enemies soi Nes To 
and ®. The resultS were compared with the value of va" calculated 
for the same input values using equation (3.10). In all cases, the 
results agreed within less than 1%. It is therefore concluded that 
in the limit of zero ionization, the Ingard~Schulz dispersion 
melations predict the correct attenuation. 
eae Neutral HWeatinusio. Electrons 
gne effect Of ehergy transteryLo the newer aie 


electrons is expressed through the term 1/0)T. Ingard and Schulz 


Setine 7 by the relation: 


1 H 
mo (YL) 5 (3 zn) 
e 
The rate of eneroay transfer per unit mass, H, is defined bv: 
ona KT 72 
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where 0, the average cross section for energy transfer, 15s 


defined as: 


: yee Fa 
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The integral over cos 9 is the total cross section 
(velocity dependent) for momentum transfer, Ap? 28 defined by 


Present Lae. ht. 1s oalven by: 


fF 


This is the correct expression for 09 __, vice the expression 
given in L7] as confirmed by private communication with M. Schulz 
(due to a typographical error). 
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a : do 


dy( Vv.) he 2 b=cos DD) eee d(cos @) (3,16) 


Therefore: 
3)" iNe 8/27? 


OnE UE I a =. v dnl )F (v,) dy, : (327 ) 
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A Maxwellian distribution of electron velocities is a 
good mak oxiniacaion to the actual distribution in a gas discharge 
only at extremely low values of E/p feo." it as recognized that 
E/p is not low in the usual case as evidenced by f19gunes 2 5eand 
126 of von Engle [40]. However the Maxwellian distribution is 
computationally convenient to use and is valid as a first approxi- 
feueron to the real distribution function. It is felt that the 
rovement in accuracy resulting from the use of a hetter 
approximation such as the Druyvesteyn distributicn does not justify 
the increased computational difficulties. The value obtained for 
us under this approximation is smaller than the actual value 
because the approximation indicates a greater number of high 
energy electrons than actually exist and the cross section de- 
creases at higher electron energies. Substitutirg the expression 


for the Maxwellian distribution into Eq. (3.17) gives: 


2 
"ert ae ay aco ae 
go =(|—< = , 
en m at eine: meer dv, 
Oo 


The cross sectron for momentum transfer iS available in 


terms of electron energy expressed in electron volts. By defining 
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ca; KT /e; the energy equivalent electron temperature, and 


effecting a change of variable, 0 on can be GxeiesceC sa nmeric 


simple form 


Lea) 
1 2 -R 
== 3 dR “AS 
en a6 chy) Sadar og ao) 
Wa, 
where R = as is the variable of integration by which the 
e 


average cross section for Sea Granster can be evaluated from 
the momentum transfer cross section data. 
5. Collisional Losses 
If we refer back ao equation (2.2), we see that the terms 
involving the constants of proportionality ee arise from the 
collision integral for momentum transfer. The collision integral 


can be written: 
N 


Chee: 
Z cine A = = be 7 
Jas, ene) 5 ) ow, De ey ve) PECs va) : 


In order to give a numerical value to W it must be expressible 


jk? 
in terms of the variables and gas properties. To do that, its 
definition must be understood. Since the momentum transferred 


from species 1 to 2 through collisions is equal and opposite to 
wae transferred from 2 to 1 through collisions, then We can 


equate; 
le. ee = age a ee 


Mars shows that the constants W.. are related to W . as: 


jk 8) 
N m 
ee wi: k 
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Notice that due to sthe mass rage tacvor. are and oe cannot both 
be collision frequencies. If in Eq. (3.19) the j-th component is 
taken to be the heavier component and the k-th component the ion 
or electron component, the resulting expressions will be con- 
sistent with the relations given by Ingard and Schulz. [In 
addition, when the heavier component is the neutral gas, this 1S 


consistent with the picture of W., being the collision frequency 


Jk 
of a particle of the j-th type moving through the k-th fluid 
component since the neutrals are much more numerous. We can there- 
fore develop analytical expressions for W , W.. andW... The 
en in el 
other terms will be related by Eq. (3.19). 
a. lon-Neutral Collision Frequency 
If it 1s assumed that the i1on-neutral collision cross 
section is essentially the same as for a neutral-neutral collision, 
then it follows that the ion-neutral collision frequency is the 
same as the neutral-neutral collision frequency. The neutral- 
neutral collision frequency 1s very closely approximated by the 


elementary kinetic theory result for a hard sphere ideal gas. The. 


1Lon-neutral collision frequency can therefore be expressed by: 


=F =\55 § 
“nin JP = uae (3-29) 


b. Electron-Neutral Collision Frequency 
in Order toOuobtain an expression for the elleetron- 
neutral collision frequency,we must start with the integral form of 
the expression given by Present [36] for the collision rate for two 


intermixed species. 
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7 a({v_-v l) 
a 3 3 ee ee f 
PEN fe Ye g Yn d Ee Veal ai) aie! 
For electron-neutral collisions, Ivo-v I= vy and) the sinveqmarlks 
can be separated. The frequency of collision of an electron with a 
Meutral particle, we is Ban No: Substituting also the expression 
for the Maxwellian distribution for fees), the collision frequency 


can be written: 





| _ mv 
N fm. 3/2 ae Ge 
Ww = fo Vv O(v Je dv 
en KT e e e 
\ 2m e O 


Finally this can be expressed in a form similar to the expression 


for 0 as: 
en 


| 2KT a ee = 

A = \ a enn Ney e ar Eth eae NAN 

rae ic \ he { VER IRE a Voce 
e O 

c. The Electron-Ion Collision Frequency 


The electron-10n collision frequency for momentum 


transfer Wis 1S given by Tanenbaum [41] as: 


~12 Tee 
. = 3.62 O j. = 
wen 2 oe Nee rey Nn (N. in cm 3) (3522) 
where A = 127 7 No is the Coulomb cutoff parameter, 
and ty? the Debye length is given by: 
I 
2 
7 KT 1 
Tr ss eee em > 
D m Ww 
e e 
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Expressions for the plasma frequencies for electrons 


me, and for 10ns W_, 
e at 


A = 
gq” N 
W = e e 
e m € 
eo ; 
oa *5 
g 
W. = a ee ww) 
Ze % e 
18 ae @ 





are well known and present no computational difficulty. 
6. TIon-Acoustic Losses 
Bosses caused by Viscosity and thermal conductivityean 
tre. 10n fluid Br oneee are represented by the term W/W). which 
-1s defined by 
irae x 


. / 
a : Fiat ve Soi) 


The quantity Via; represents viscous stress in the ion fluid. The 








Ele 


oiner term Fos. when multiplied by e.; the ion temperature per- 


feation, accounts for heat conduction in the ion fluid. Adapting 
a ae Approach, assume that w/t), has the same form as 

W/W Given by equation (3.9). Under such an assumption, it can 

be expressed as: 


L 


= [2.4 w ra +e 


Per ee 
; \2.4 Wh .c. (3.23) 





c 
Loie 


sgl 
where the subscripts indicate that these are ion terms. 
It 1S not unreasonable to assume that, since their masses 


are essentially the same and their temperatures nearly the same, the 


49 








ion sound speed Cc. 1s the same as the neutral sound speed C. 
However the mean free path for ions Le, must be unlike that for 
neutrals, Since one must consider the effect of Coulomb inter- 
actions. wu [42] expresses the mean free path for Coulomb 


interactions, which have been corrected for Debye screening, as: 


1 


LL. = eonw__— F (ise) 
4 ON or” &(r_/r_} 
e m ‘PP m 
2 
= e e 
where Une Kt 1s the mean distance of closest approach for 


two like charges q. Wu refers to T as the mean plasma temperature. 
Since we are considering i0n losses, the temperature to be used in 
the expression for x is the ion temperature. However the Debye 
length is determined by the electrons due to their higher tem- 
perature and greater mobility. 

Because the mean free path of ions in a weakly icnized 
gas is limited not by Coulomb interations but through collisions 
with neutrals, this expression alone would give erroneous results. 
The mean free path, 4, can be expressed in terms of the mean 


velocity v and collision frequency v as 
pe ye 
ae 
In a mixture of gaseous components, the total collision 


frequency of one component is the sum of the partial collision 








frequencies. This leads to the result: 
tg ee 
\ oft ie Nad May, 
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ea .. and making the approximation thatA. = the 
Ignoring 1A . nd making the app m n Ae eS 
mean free path for ions can be written as: 


2 1 = 
ee i: We aye : | (3-25) 


B. COMPUTATION AND EVALUATION 
len “Defining the Limits of the Variables 

Since the goal of the analysis is a prediction of results 
for a given experimental arrangement, it would have been desirable 
to formulate the problem in terms of filling pressure Po? iCall lag eve 
temperature aoee tube radius r, frequency v and current-I.. Unfor - 
tunately the conversion of Po (eee and Ty to the variables To TY 
and ie is not a trivial problem. For that reason, the independent 
variabies chosen were: neutral density Nee neutral temperature Vy? 
electron density Nee electron temperature To and angular frequency 
wW. The limits of these variables were then chosen to correspond 
to the values expected in a feasible experiment. 

A practical range of values for the tube radius is from 
-5 cm to 10 cm. This value-was fixed at 2 cm except when investi- 
gating its effects. 

The lower limit of the Sees frequency was fixed at 
100 radians/second because it was experimentally possible and 
numerically convenient. The assumption of plane waves was implicit 
in the formulation of the problem by Ingard and Schulz. From an 
experimental point of view, in order to avoid the effect of higher 


modes propagating in a tube, this sets an upper limit for W of the 


sul 








order of Done ie ih 25) - However most pulse experiments involve 
higher frequencies and yatsul, et.al. Ua worked with 50, 100 
and 200 kilocycle continuous and burst sineisoidal signals. The 
limit of validity of the continuum description is when the wave- 
length is of the order of the mean free path. This corresponds to 
about 10? radians/sec. at room temperature and pressures around 

1 torr and was selected as the upper limit for W. 


ewe 


fhe neutraledensity Pimitsy Were Ser at 325) <0 oemem 


18 -3 


fee. > XxX LO cm ” , corresponding approximately to filling pres- 


sures of .Ol torr and 50 torr at room temperature respectively. 
The electron-density limits relative to the neutral 


density were selected arbitrarily, with a logical basis, as 


197 on to 107 3x with absclute limits, of 10 a 


nN nN 


To aoe 


£00 15s 10 am 

For currents of milliamperes at low pressures around 1 torr, the 
degree of lonization is about ome [40,43]. aA lower limit two 
orders of magnitude below this is reasonable and a higher limit of 
three orders of magnitude greater would correspond to ampere cur- 
rents. These relative limits were tempered by the current densities 
which they implied for large pressure changes. 

The neutral gas temperature is not truly independent 
but actually depends on the discharge current, pressure, tube 
EBeadius, and external cooling effects. For the purpose of this 
analysis, it is sufficient to assume the temperature can range from 
the boiling point of the gas to the upper limit of temperature 


which can reasonably be obtained in an oven, approximately 700 Kelvin. 


ae 








The electron temperature Te 1s definitely not an in- 
dependent variable, being determined by the discharge conditions 
specified through the other variables [40]. The only analytical 
expression available implies that the electron temperature is a 
function only of the gas and the similarity parameter pR. The 
expression is of limited validity, falling at very low values of 
meeand high p. Rather than restrict Te with an artificial con- 
straint, it was felt that treating it aS an independent variable 
would yield more valid results. The limits were chosen to cor- 
respond to those normally encountered experimentally [40,44], 


4 To 10° 


10 Kelvin. 

It should be pointed out that permitting these 
variables to move freely within their limits places less con- 
straints on their possible values than actually exist. The re- 
sult may predict physically unrealizable results. 

2. Computation of Aitenuation and Determination of Minimum 
For the computation, the formulas were first placed in the 
order in which their values were needed. They were: then translated 
for computer processing and a program written using the FORTRAN IV 
programming language. The values for the fundamental constants 

are from LAs), data for Argon are from Cook [46] except the atomic 
diameter and scattering cross section for momentum transfer which 
are from McDaniel [39]. 

The scattering cross section data were taken from figure 


4-3-1(c).of McDaniel, using a transparent overlay grid and dividers 


me establish a table of values of a 


iy Versus electron enera in 


elvcetron Volts. <A subroutine was written to use this table to 


og 








perform the numerical integration to determine Ose using Eq. (3.18) 
and to determine Oe using Eq. (3.2) jp. in) Compniging ae the 
momentum transfer cross section data were used to approximate the 
collision cross section. The momentum transfer cross section is 
from 20% to 50% smaller than the collision cross section for all 
values of electron energy. This substitution causes the electron 
collision frequency to be low by ~ 25%. If this was found to be 
critical, a correction could be applied. This procedure was con- 
sidered to be justified by the increased convenience in setting up 
the problem and in reducing the number of complete integrations 
required by @ factor of two. 

16 10 


NaS eee 


A set of fixed values: W = 10°, N = 10 3 
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of toeevalwate all expresoquems wath 
an electronic calculator. This served as a cross check on the 
computer results. To check the validity of the expressions used to 
calculate attenuation, a comparison was made of the attenuation 
factor for bulk losses with a value from another source. Using 


baye (3.9) and (3.7), we can write 





a 
BULK _ 1 L Z a 
B= os (3 = |G es 
Ze (@: 


n / BULK 


therefore 


“>a 2 £ 
a aes: 


Evaluated at standard temperature and pressure, this yields 
q = ~ 
ae 1.89 x 10 oe 

fe 


mean 








: Z - 
mhach agrees wath the valwenot.cc) «6 io 7 sole) ae for Arco 


mm Kinsler and Frey [47]. 
32. Discussion of Results of Calculations 

At the point used for checkout of the computer program, 
the agreement between the attenuation constant calculated using 
Eq. (3.1) and the Ingard-Schulz approximation, Eq. (2.21), was 
very good with differences on the order of a few parts per million. 
However the agreement was not so good in other areas. The use of 
Eqn. (2.21) yielded predictions of amplification over a wide range 
SEeconditions. In contrast, results using equation (3.1) show no 
amplification sine the entire range of the variables, but merely 
an increase or decrease in attenuation, with respect to that with- 
Suite oOnization. Phe amplification predicted by Eq. (2.21) a 
Mienerore considered erroneous, arising from factors being neglected 
in the approximation. 

The minimum attenuation was found to occur at the limits of 
the allowed ranges of the variables. Those variables at their high 
limits were neutral density, electron density and electron tempera- 
foube. The FeGueney sane neutral temperature were at their low 
iamits . 

Figure 1 shows the attenuation ae AatMnCctE Von Of anaulanr 
frequency, displayed on a semi-log plot. The curve parameter is 
neutral density. Neutral temperature, electron temperature and 
degree of ionization have been held constant at their optimum values. 


A more revealing display is shown in Figure 2 which is a plot of the 
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Figure 1. Predicted Attenuation Factor in Discharge 


Plasma in Dependence on W for Various Values 
of Neutral Density. Parameters are: fee S&K, 
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Relative Change in Attenuation due to Ionization 
in Dependence on W for Various Values of Neutral 
Density. 


Values calculated using Ingard-Schulz equations 
with parameters ee 88K, La 109K, and 

ee = ioe The values of neutral density 
for the curves are: 


1. ve =” Tees i018 ae? 


2 i 10.7 ene 
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By 








fractional change in attenuation due to ionization. The curve 
parameters are the same as for Figure (1). Positive values cor- 
respond to an increase in attenuation; negative values, a decrease. 

The ion perturbation velocity was found to remain the same 
magnitude as the neutrals and in phase. The electron perturbation 
velocity was the same magnitude as the neutrals, and in phase, at 
low frequencies, and decreased to zero and went out of phase at 
high frequencies. The frequency at which the transition began 
shifted to higher frequencies as the density increased. This is 
displayed in Figure 3. 

The results of the computation indicate that net amplifi-~ 
cation is a false prediction, arising from approximations; however, 
attenuation is aftected to a measurable degree. the inaximum dGe- 
crease occurs with high pressure, low temperature, low frequency 


and a high electron temperature, highly ionized discharge. 
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IV. EXPERIMENTAL AKRANGEMENT AND PROCEDURE 


The original concept was to use an auxiliary discharge, with 
a gated sine wave of current to produce gated sine waves of sound. 
the sound amplitude was to be measured before and after passing 
through the main discharge area. From a comparison of the signal 
amplitudes at the two points for various frequencies, the attenu- 
ation factor and its variation with current magnitude and polarity 
would be determined. The analysis of a preliminary experiment 
showed that other factors had to be considered and that a Fourier 
analysis of the Signals was needed. Compromising between the need 
for a signal with low frequency components and one with all signi- 
ficant features appearing before the reflection of the signal from 
the tube ends, a one millisecond current pulse was selected to 
provide the sound pulse. Although a two ampere power supply was 
available fox the main discharge, the high temperature of the main 


discharge column limited the experiment to half an ampere. 


Pee SOUND SOURCE 

The sound source was a smali auxiliary discharge far removed 
from the main discharge and the two diodes bounding the main dis- 
charge. The auxiliary discharge was always on when in operation. 
The quiescent state was a normal glow with a current of less than 
10 ma. The peak current during a pulse was around 150 to 200 ma. 
which corresponded to an abnormal glow. The current in the auxiliary 
discharge was controlled through a transistor load. During pre- 


liminary work, it was ascertained that the sound produced «as 
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elosely related to the current in the auxiliary d%eecharces utenor 
to the voltage or to the power dissipated. Therefore only the 
current was controlled and monitored. 

Although a pulse was used for the final data, the current con- 
trol circuitry was designed for sinusoidal currents, either con- 
tinuous or gated. To avoid extinguishing the discharge, a minimum 
current flow was required at all times. This required a bias signal 
equal to the sinusoidal signal. The continuous mode was not opera- 
tional due to rapid deterioration (in times less than or of the 
order of one-second) of the discharge current hen run at the high 
currents required for a good acoustic Signal. The aes was not 
thoroughly investigated, but was believed due to excessive power 
G@iestrpation. 

The current control circuitry is shewn (enclosed in dashed 
Hines) in Figure 4. It consists essentially of two circuits feeding 
the base side of two transistors operating in parallel. One cir- 
cuit consists of a small d.c. battery with a variable resistor for 
controlling the quiescent state of the discharge. The other cir- 
cuit sums two input Lee eG PrOVide the a,c. sigpal, 

Because of their high voltage capability (300 V), and their 
high dynamic resistance, Motorola MJ3202 transistors were used. 
Parallel operation was chosen to give a greater current range to 
the discharge while remaining within the power rating and linear 
range of the transistors. Emitter biasing resistors were used to 


balance the load between the two transistors. Although some 


Gl 












GRIEZI eG 
PUGCE 
GENERATOR 


GR 1396A 
TONE BURST 
GENERATOR 






Figure 4. 




















300 VOLT 

REGULATED 
D.C. POWER 
SUPPLY 





Sound Source Circuit. 


62 












negative feedback occurred due be Ne absence of bypass capacitors, 
it did not have a noticeable effect. | 

A Generai Radio GR1217-C Unit Pulse Generator provided the d.c. 
pedestal for the gated sine wave. The degree of modulation by the 
gated sine wave and the amplitude of the a.c. discharge currents 
were adjustable through variable resistors. For ease of control, 
multiturn potentiometers were used for control of bias, quiescent 
Sarrent and “2c. Current amplitude. 

The General Radio equipment used required the low side of the 
transistor circuit to operate at ground potential. To minimize 
electrostatic interaction with the diode microphones, the anode of 
the auxiliary discharge was grounded. This meant that both ter- 
minals of the d.c. power supply were floating above around. Since 
the floating potential varies greatly during the current modulation, 
capacitive coupling of the terminals to ground was found to distort 
the discharge current waveform. The power supply used (Universal 
Electronics Model G3050M) was chosen because of its low ccupling. 
However its maximum voltage was 300 volts. This did not permit 
using the full current capability of the circuit at low pressures. 
Consequently operation was limited to pressures above 5 torr for 
good acoustic signals, although use of another power supply would 
have extended the operating range. Figure 5(a) shows a typical 
discharge current pulse waveform with a gated sine wave. 

The electrodes for the auxiliary discharge were two tantalum 
Bybinders (17; cm long by 2 cm dia.) with faces Roar by 35 em. 
They were mounted axially in the tube, supported by .070 inch 


diameter KOVAR leads. The leads were beaded with glass for electr:- 


el insulation. 6 3 








(a) 


(b) 





Figure 5. 


Sound source current waveform and resultant 
sound. Ambient pressure is 8.7 torr. Gas 

is Argon. Frequency of sine wave being gated 
is 3000 Hz. Time scale is .5 ms/cm. 


(a) Sound Source current. Vertical scale is 
100 ma/cm. 


(b) Sound produced by discharge current. 
Vertical scale is uncalibrated. The second 
pulse of sound is the reflection from the 
near end of the tube. 
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Be SOUND DETECTION 

Mobility-limited thermionic diodes were used to detect the 
sound waves. They were mounted axially, one on each side of the 
discharge region. Electrostatic shielding was provided by a 
cylinder of heavy, 2.3 mil., copper f011 which was wrapped around 
the tube at the location of the diode. The shield extended several 
inches on both sides of the diode and was connected to ground. 
ince the cables and connectors were shielded, the only part of the 
circuit not shielded was a small portion of feedthroughs where 
they penetrated the glass. Also the diode was not shielded on the 
ends. For a more complete discussion of the diode and its sus- 
ceptibility to noise, see Appendix A. Figure 5(b) shows a typical 
diode response signal. 

For the experiment, filament power for both diodes was provided 
by the same power supply, a Power Designs 36100R regulated power 
supply operating in the constant voltage mode. Figure 6 shows the 


circuit employed for power and signal detection. 


See lUBE DESIGN AND CONSTRUCTION 

There are two features of the tube used in the experiment which 
are worthy of note. One is the method of mounting the elements. 
The other is the length of the tube and the placement of the elements 
therein. 

The tube elements were mounted on feedthroughs which did not 


distort the tube walls significantly. The feedthroughs were 
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sealed with a special glassblowing technique which is discussed 
in the appendix. The result was a tube which could be treated as 
an acoustic waveguide for frequencies much higher than those used. 

To reduce the effects of reflections and electrestattemmnter— 
actions, the elements were physically separated so that the acoustic 
time of flight separated the unwanted, from the wanted, signals. 
This required a half meter buffer region at the tube ends and a one 
meter region between the sound source and the nearest diode. 

In addition, it was discovered during preliminary experiments, 
that there were positively and negatively charged particles ex- 
tending into the supposedly neutral region far beyond where they 
were expected. A series of experiments established that they were 
real and that hoth tyvnes were present. The buffer region reduced 
the effect but screening electrodes were used for further reduction. 
eanee with cold filaments, any diode signal — attributable to 
charged particles rather than sound, the screnning effectiveness 
was studied with cold filaments, in which case, the diodes acted 
like double Langmuir probes. The cold filament effect was eae 
in Neon but barely noticeable in Argon. Although interesting in 
itself because these signals were detected at distances from the 
discharge of 30, 40 and 50 cm in Neon at pressures ~ 10 torr, the 
phenomenon was studied only encugh to eliminate its interference. 

A schematic drawing of the tube appears in Figure 7 and the lo- 
cations of the tube elements appear in Table 1. 
Copper-cconstantan eReUnocoup lee of 5 mil wire, were used to 


measure the temperature on the tube axis at three points. For 








TCI TC2 T.C3 


Al = A2 Di Si C1 C2 $2 D2 


Figure 7. Schematic Drawing of Discharge Tube. 
Identification of Symbols and Locations 
are Shown in Table I. 


TABLE I 


LOCATION Or TUBE ELERNENTS 


ITEM SYMBOL LOCATTION*¥ (cm) 
End ] Tt ~-149.7 
Anode i Al LOZ ah 
Anode 2 A2 -98.7 
Thermocouple 1 TEL 50S 
Diode ] Di O 
Screen 1 Sl 2053 
@irnode/Anode 1 eu) 37.3 
Thermocouple 2 EZ 102 
Cathode/Anode 2 EZ 122.2 
Screen 2 S2 oes 
Thermocouple 3 Fe3 Eo eS, 
Daocde 2 D2 L528 5 
End 2 203 + 


- Location given is location of center of element 


7 Tube ends taper from 4 cm to 1 cm diameter over 5 cm. 
Distance refers to beginning of taper. 


+ For calculating tube volume, the effective length is 363 cm. 
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Support and insulation, the leads were first sealed in fine twin 
capillary glass tubing with ULTEK VACSEAL epoxy. The wire to 
extend into the tube was left exposed. During construction of 

the main tube, a small piece of tubing, large enough to hold the 
capillary tubing, was attached. The capillary tube was sealed 
inside of this tube with the epoxy. The seal was good and the 
thermocouples worked well. Such a method 1s not recommended for 
the main discharge area. This 1S because the temperatures en- 
countered in the vicinity of the discharge may exceed the maximum 
recommended for the epoxy, 150°c. In this experiment, such heating 
caused contamination of the discharge and limited operation of the 
main discharge to currents below .5 ampere. 

The pumping train used for evacuating the tube consisted or a 
VEECO Model EP2Al1 diffusion pump, rated at 85 £/sec, pumping to a 
DUOSEAL forepump. A NRC model 720 ionization gauge control with a 
type 536-SK ionization gauge was used to monitor the pressure above 
the trap, which was kept filled with liquid nitrogen. A one inch 
line was used between the trap and the high vacuum isolation valve 
(a GRANVILLE PHILLIPS type C value) used to connect the discharge 
tube to the pumping system. The ultimate vacuum obtained was 
around 5 x Wor torr as measured by the gauge. 

The gas filling system was connected to the tube via a (CEC 
type 25) high vacuum valve. Pressure was measured with a manometer 
filled with Octoil S. A trap located between the filling system 
and the tube was partially filled with liquid nitrogen whenever the 


connection to the tube was open. 
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D. MAIN DISCHARGE 

In order to determine the effect of polarity = theserceeueac. 
at both ends of fhe main discharge needed to serve as cathode or 
anode. Therefore they were of identical design, consisting of a 
filament and a cylindrical anode. The filament was used as the 
motecathode for the main discharge or as the hot cathode for a 
"Pupps discharge" at the anode of the main discharge. It was a 
faeioch diameter, 7 turn helix, directly heated tungsten dispenser 
cathode designed for laser work and produced by SPECTROMAT. It is 
rated at an emission current of Baar at 50 coe The anode 
cylindrical electrode which was mounted concentric with the fila- 
ment, was located so that the first two turns of the filament, 
on the discharae side, were exnosed. This was done so that the 
Pupps discharge would tend to go toward the leading edge of the 
anode, and thus eliminate anode CGE Sere A small low 
Yipple (20V, 1.5A) D.C. power supply was used for the Pupp's dis- 
charge power source. 

Filament power for the main discharge cathode was Si aen eas 
a KEPCO KS-8 100M power supply. Filament power for the Pupps dis- 
charge cathode was supplied by a POWER DESIGNS PD36100R unit. The 
load for the main discharge was two high power variable resistors 
connected in parallel. The main discharge power supply was a 
KEPCO 770B, 600 V d.c. supply. 

A reasonably simple, yet effective system was used to elec- 


trically isolate the main discharge system from the remainder of 


mete Circuitry. The 60 cycle input power for all units was isolate 
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using transformers. All equipment for the main discharge was 
mounted on the same rack, which contained no other equipment. The 
chassis ground for all of this equipment was connected to a common 
point which was then connected to ground via a one megohm resistor 
Meparallel] with a 10 ff capacitor. The a.c. and d.c. potentials 
across these points were monitored with high impedance devices. 
This precaution was taken to prevent the buildup of a high elec- 
trostatic potential between the two ground points. The a.c. current 
in the discharge was monitored by using a current probe with a 
sensitivity of 1 mv/ma. To prevent a large potential difference 
between the two ground systems, arising from interaction of the 
fields, a potentiometer was used to vary the potential of the main 
electrodes with respect to the secondary ground. This is sown 
in Figure 8 as part of the main discharge circuit. 

The discharge was ignited with a Tesla coil in conjunction 
with a copper wire wrapped around the discharge tube between the 
anode and the cathode. This was necessary because of the high 


breakdown voltages required at the pressures used. 


E. ELECTRONIC MEASUREMENTS 

The three waveforms of interest - the sound source current, 
and the two diode signals - together with a trigger signal were 
recorded on a four channel FM tape recorder (Precision Instrument 
SO. , PI-6200) for later analysis. Recording was done at 37.5 ips. 
Record and playback specifications are contained in Appendix C. 


The frequency response did not appear to distort the two diode 
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Signals but it did affect the trigger and discharge current signals 
since these had significant high frequency components. Figure 9(a) 
shows a typical diode signal. Figure 9;:b) shows a typical playback 
of a diode signal. It should be noted that these figures are not 
derived from the same sample. 

The trigger signal was provided by a General Radio GR 1217-C 
Unit Pulse Generator operating at a pulse repetition rate of 
10 Hz. The output signal was a one millisecond rectangular pulse 
which was used to trigger the scope and was also recorded. The 
"delayed sync" output was used to trigger another unit pulse 
generator from whies a one millisecond pulse was fed to the cur- 
rent control circuit. The discharge current was observed by 
measuring the voitage across a 5 vim precision resistur. 

The output signals from the Giese were amplified with two 
differential a antite ee TEKTRONIX RM122, with an input impedance 
of 10 megohms. The amplifiers could provide amplification of 
either 100X or 1000X. Both output signals were recorded on tape. 
The diode output contained very sharp noise spikes which, after 
Bitication, adversely affected the recorder operation. To 
eliminate this, a slicing circuit, which clipped all signals whose 
Pemue exceeded + 1.5 volts, was employed between the amplifier 
output and the tape recorder. This circuit permitted an improvement 
in the signal to noise ratio since the noise could not exceed 
1.5 volts and the signal could be amplified to a level just below 


that. 
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Figure 9. Comparison of Typical Diode 
Signal with Playback of 
Diode Signal. 


(a) Direct display of diode signal. 


(b) Playback of diode signal. Due 
to triggering problems, the 
playback is not a playback of 
(a). The playback does show 
how little the recording process 
distorted the signal. 
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F. EXPERIMENTAL PROCEDURE 

After preliminary experimental work had established the back- 
ground, the progress of the experiment can be divided into five 
Buages: (1) construction of the tube, (2) preparation of (thestube, 
(3) calibration of the diodes (4) data accumulation, and the last 
stage, (5) analysis of results. The analysis of results will be 
discussed in the next chapter. 

ie GOnNStruction 

Due to its great length, the tube was constructed in three 

sections. After construction and before annealing, each section 
was leak tested with a helium mass-spectrometer. The sections 
were annealed in a vacuum and then sealed under vacuum until time 
for assembly. The supnort svstem built for the tube consisted of 
a long aluminum beam with six rigid tube supports. This permitted 
Pemsections to be rotated while being audi together so that the 
joint was reasonably smooth with no sags or large irregularities. 
After the sections were joined together, the tips of the thermo- 
couple support tubes were broken, the thermocouples inserted ae 
sealed with epoxy. The epoxy was then cured for two hours with 
heat lamps. After curing, the tube was placed in its final 
position and sealed to the system via a section of glass tubing 
containing a metal bellows. 

Zor epeatation 

Preparation of the tube involved bakeout and outgassing 

of the tube and its aienente so that a high vacuum could be ob- 


tanied followed by filament activation. Bakeaout was perfc ‘ied 


Te 








with heating tapes at 200°C for over 24 hours with the system 

under evacuation. The cold trap was kept filled with liquid 
nitrogen following bakeout. Outgassing was performed by inducticn 
heating, after bakeout. Since the leads were KOVAR, with a much 
lower melting point than the tantalum cylinders, outgassing of the 
tantalum at elevated temperatures was done in short bursts. In 

the case of the diodes and main discharge electrodes, the filaments 
were activated first and then the anodes were outgassed. Filament 
activation consisted of slowly raising the temperature by varying 
the filament voltage until the temperature, as measured with an 
optical pyrometer, was at a certain value. For the diode filaments, 
this temperature was WO Coe (subscript B denotes brightness 
temperaturejtine .:ecommeided temperature for thoriated filaments. 
Following the manufacturers instructions, the main cathode fila-~ 
ments were Only heated to UIBON Ce This temperature was held for 
five minutes. 

Following outgassing and filament activation, the filaments 
were aged for several hours at their normal operating temperatures. 
Then the main discharge anodes were cleaned by electron bombard- 
ment. This was accomplished by applying 300 volts to the anode, at 
which time the anodes turned bright red under the absorption of 
45 watts of power. 

Be Calibration of the Dicdes 

After the tube was prepared, with all elements outgassed, 

the diodes had to be calibrated. Since the mobility-limited cur - 


rent is less than the space-charge limited current, and th« 
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transition is a function of gas pressure. it is best first to 
attain space charge limited operation. This was easily determined 
by preparation of a log-log plot of the plate current versus 

plate voltage to a value of plate voltage (30 volts), well in 
excess of that to be used (13% volts). This was found to be a 
straight line, indicating that the filament voltage selected 
provided a good operating point. 

Calibration consisted of varying the pressure in the tube 
and measuring the plate current as the plate voltage was varied. 
The filament voltage was held constant. As the pressure increases, 
there is increased filament cooling. It is not practical to use 
an optical pyrometer to maintain constant filament temperature at 
all pressures. Since the power applied increases as the tempera- 
ture decreases for constant voltage, the constant voltage mode is 
semli-compensatory. 

The filaments qeue thoriated tungsten and at first operated 
extremely well. After the main discharge was first run, which 
Overheated the epoxy seal at the thermocouple, the diode current 
became noisy and started to decrease. [t became necessary to 


increase the current to a high value and operate the diode fila- 


ments aS pure tungsten. This required a recalibration. 
4. Data Accumulation 
It was decided to take data at three pressures. The low 


pressure was selected to be around 8 torr because a larger power 


supply was needed to obtain a good current pulse at pressures much 
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lower than that. The high pressure was 16 torr because of the 
marficulty in getting the main discharge started at higher 
pressures. The third pressure chosen was 12 torr. 

In order to be able to detect the existence of any trends 
(in the data), due to time the order of experiments was chosen 
with a random number table. First the order of the pressures was 
chosen and then the order of the polarity of the main discharge. 
ihe order of runs is shown in Table II. Due to the difficulty 
of filling to the same pressure, the runs were not completely in- 


dependently crdered. 
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A set procedure was followed when taking the final data. 
First all equipment was turned on and allowed to warm up. Then 
the tape recorder was calibrated. The calibration of the oscil- 
loscope was checked. For each pressure, the steps were the same 
except for minor variations due to the different order of the 
polarities within each pressure. The trap in the filling line 


was partially filled with liquid nitrogen. After the tube was 


filled with Argon at the desired pressure, the valve to the filling 
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manifold was secured. The diode and main discharge filaments 

were given plenty of time to warm up. When the main discharge 

was on, the Pupps discharge was adjusted for minimum AC current 
variation in the main discharge. This was observed on the 
oscilloscope using the current probe. The d.c. diode signal was 
measured after adjustment so that the amplified a.c. signal did 

not exceed the clipping level of the slicing circuit. When 
everything was ready, the tape recorder was started. While it 

was running, the temperatures were recorded. Between runs, the 
tape recorder calibration was checked and the output zero adjusted, 
but the input gain was not touched. After the experiment was over, 
a 5000 Hz tone burst sine wave, stepped down to 2 millivolts peak 
tO peak was fed to the amplifiers im parallel. This was recorded 
for later use as a calibration aid. After each set of runs at a 


single pressure, the tube was pumped down to a hard vacuum and 


meal tecd for the next set of runs. 
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V.. THEORETICAL TANMALY SiS 


This. section treats the analysis of data theore treaty ai 
develops the expressions needed to utilize the data to obtain 
information from it concerning the processes under investigation. 


In doing so, it takes into consideration the progagation of errors. 


A. ENSEMBLE AVERAGING 

Nine experimental runs were made, each with different conditions. 
An event, or experiment, consists of a current pulse, and the sub- 
sequent sound production. An experimental run consists of many 
events or replications of the experiment. The data recorded for 
each event were the triggering signal ard the waveforms at three 
points. After conversion tu diyitai form, each waveform was te- 
presented by a set of samples. Since a run consists of a set of 
events, the data consists of a set of sets of samples which shall 
be referred to as an ensemble of sample sets. 

Consider just one of the signals. The samples in any one of the 
sets were not taken randomly , but have faxed relationships. cto each 
other which reflects their time dependence. For the set of samples, 
the zero of time is chosen to correspond to the triggering event. 
Thus, the n-th sample represents the value of the signal at time 
t = nT where T is the period of the sampling frequency. 

This sample represents not only the time signal at that time, 
baue Gul something else. It represents the value of the cross 


correlation between the signal at that time and a Dirac delta 


80 








function signal which coune ida with the triggering event. This 
is a subtle point which was used for, advantage in eliminating 
noise. 

Mie Cross correlation function /(x) of two real funetiens 


g(x) and h(x) 1s given by: 
f(x) = g(u-x)h(u)du 


eteqg be the Dirac delta function 6(t) and h be the signal s(t). 


( 6(u-t)s(u)du 


Then 


J(t) 


and therefore, 


J(t) = s(t). 

By ensembie averaging the sets of sampies, we ovbiain tie 
ensemble average of f(t), which eae cross correlation function 
of the signal with the-triggering event. For a large enough 
ensemble, this should approach the time average of the cross cor- 
relation function. In the time average, contributions from non- 
@errelatable portions of the signal, such as those due to noise, 
random or A eeCony will vanish. 

Because of their intimate relation, f(t) and s(t) may be used 
interchangeably except that it must be Bayamaniizied that in cach 
sample set, s(t) has a different reference time while f(t) has the 


same reference (for this reason, s(t) will be used for f(t) except 


when the speciai features of f(t) are needed). 
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When f(n7T) is ensemble averaged, the mean value and its 
standard deviation are obtained. Since there is no reason to 
assume that the distribution of noise or fluctuations in value 
of f(nT) should depend on n, the standard deviation is considered 
a constant. Then, if i (n7) represents the value for the i-th 


set of samples, the process cf ensemble averaging can be repre- 


sented by 
1 N 
fry He i oa) 
( N I=- 1 
(5.1) 
N 
Z 1 2 Z 
Le) = ——= ; T)- T 
n N=2 ae Ss (a ae ) 
and 
M 
Je ys (5.2) 
M et Lh nN 
where M is the number of points in the sample. It is possible 


to consider this ensemble average f(nT) and the standard deviation 
@ as the expected value of the signal, s(n7) and its standard 


deviation, O. 


B. FOURIER TRANSFORMAT IONS 

When a harmonic waveform is fed into a linear time-invariant 
device, the output is also harmonic but in general with a dif- 
ferent amplitude and phase [48]. This is true whether the device 
be a filter, or transducer or some other system with an input and 
an output. Such a device can be completely specified by a fre- 


quency dependent complex quantity Z(vV) given by the ratio of the 
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eitput to the input. Sit sthewinpwet v(t) is given by 


WALA» 
Vea ee Ae~ 
li 
and the output by 
L20VT+i 
eee ena en 
2 
then 
B if 
=o ee 
A) a 


If an arbitrary input v,(t) 1S applied to the device and it 
is desired to calculate the output vo(t), Uae) is analyzed into 
its spectrum; each spectral component multiplied by the corre- 
sponding transfer tactemmto Obtailmene spectrum of the output, and 
then v5(t) is synthesized from its spectrum. Thus 


V.(¥) = 2(v)v,(v) 


and 


Poy t 


v,(t) = ie Z2(v)V,(v)e dv 


(Here and elsewhere, where a quantity is viewed in both the 

frequency and time domains, the upper case symbol will be used 

for the frequency domain and its lower case for the time domain.) 
If the process of convolving g(x) and h(x), is represented by 


g*h, defined by a 
J = 97h ={ g(u-x)h(u)du, 
(oo) 


and the process of Fourier transformation is represented by 2 , 


the determination of v(t) can be presented by the diagram below. 
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v,(t) 3 Vy(V) 
* x (mul tip 2a cations 
Zt) Si) 


l, 


va(t) > V3(v) 


In this experiment, there are three time dependent waveforms, 
mec). s,(t) and S(t), corresponding to the sound-producing dis- 
eierge current, the sound at point 1 and the sound at point 2 
respectively. The object of the experiment was to obtain the 
Geaaster function of the pabuer saan Cf CUrPenErto Sound ZA(V)s 
and the transfer function of the spatiai region between points 1 
and 2 Zs5(V)- The entire process can be represented by the 


diagram below, where the time representation of the transfer 


functions and the convolution process has been omitted for clarity. 


bc IT(v) 
x 


z,() 


.¥ 


S4(t) So(v) 
Oe 


Zoi () 


Sy (¥) 


V 


s,(t) 
x 


Z19(V) 
| 


Ww 


Sj(t) > S$,(v) 
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The application of Fourier transform theory to this problem 
assumes linearity of the processes. That is, if S5,(t) is the 
Setout at point 2 tO an input Of 1,(t) and if S55(t) 1s the 
response to in(t);, then the response to i(t) = 1,(t)+1,(*) is 

7 fs 7 
S(t) $5,(*) S5o(t)> regardless of the choice of 14 and 1, 


In the transform domain, the process from current to sound at 


point 2 is represented by the set of equations: 


So(v) = I(v) Ze(v) 
S)(¥) = So(v) Zp4(¥) 


S,(v) = S,(V) Z45(¥) 


Since i(t), S,(t) and S5(t) are known, their respective 
transforms are known. Therefore these constitute tnrvec equaticne 
with four unknowns, Sol); Zio(V); Zo1(¥) and Z.(V)- For the case 
where the main discharge 1s on, only ZO 1s desired. This 


can be solved for easily, 


[Sy 


To determine Ze» when the discharge is off, acoustic theory 
provides a model from which a fourth equation can be obtained. For 


sound propagating ina tube, the transfer function between two 


es : 
meimtS 1S given by ik(v)d,, 


Z15(V) =e 


where d,5 1s the distance between 1 and 2. With this equation and 


the additional assumption that, in the absence of a discharge in 
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the intervening space, k(vV) is the same throughout the tube, the 


system reduces to two equations in two unknowns, 


ik(Vv)x, 
I(v) Z.(v)e 


Si (v) 


ik(V)x,, 
S5(V) sa DKS) ) Z.(v)e 
The propagation constant k(v) is complex, 


k(v) = @(v)+ 1B(y). 


Under the experimental conditions, boundary losses dominate. In 
times case, B(v) = BOY - If phase information is neglected, 
meaning that the delay between 1(t), S,(t); and S,(t) 1s neglected, 


then By can be obtained from Si(Y); S a(V) by 


NK 


a 
A 5 SAE 
em ae 
See eT 
Iso(v)| 
Ir(v)} 


and 


I 


Izy) 


See ACOUSTIC TRANSFER FUNCTION 
After ensemble averaging, the two diode signals are represented 


by sample sets and standard deviations 


$,(nT), o 


and 
n S,(nT), Oo. 


Since the standard deviations are considered time independent, each 


BepresentS a constant signal in time, the transform of which, 2, 
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is a zero frequency component of the same magnitude as 0. There- 


fore, in the transform domain, the two diode spectra are 


Svs) 
and So(v,) 
Zz 2 
where ay and 2, affect only S(v_ = 0). 


The standard deviation or error in S,(Y;3) and So(v5) is that due 


to the noise, plus that due to other sources. (The error in the 


value of a quantity Such as S, is represented by A preceeding the 


iL 


Eeimeel for the primary quantity, e.g. AS) +) 


Let Y(v;) represent the natural logorithm of Zi5(V;)- 





Then FSA): 
MAE i ae (5-3) 
PVN’ 579 
and ‘ : 
AS, . AS, 2) % 
= — = + - 3 4 
AY(V,) é S ea) 
2 1 


For the discrete case, the earlier expression for By 1s not very 


useful. Another expression which can be written containing Be 1s 
2S ; 
Y(V5) == Poly )Yy (5.5) 
This expression can be solved for Pat each value of V.. However 


a better approach is to use all the information available and 
minimize the error in Bo: If all values of Y(v;) were known to 
the same degree of accuracy, then a simple linear least squares 
fit would provide the optimum estimate for BS: However, from 


the ensemble average and another consideration to be discussed 
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later, the error in Y(v5) is known. A weighted average will 
therefore improve the estimate for B.°- The values of Y(v5) are 
best weighted with a weighting function 


Ww, = (5.6) 


Zz os 
AY(V. ) 
where A iS a normalizing constant given by 


A = a AY(\; ) 


il 
D. SOURCE TRANSFER FUNCTION 

When there is no discharge between points 1 and 2, it is 
assumed that the transfer function is the same throughout the tube. 
This means that the value of es may be used to extrapolate the 
sound Spectrum observed ai points 1 and 2 to the origiti. There it 
may be compared with the current spectrum of the sound source 
@alcemarge to obtain the source transfer function. 


The spectral components of Sj and S, projected to the origin 


are given by: 


al 
@ 


He a8 
S5.(V4) = S;(v)e a (S.7) 


ee 
(j=1,2) 


If BO PwUsed, the impiveagion 1s that, the estimate of B is 


independent of the 2a But since B was derived from the estimates 
Or S5 and Sy this is not true. Recall that 
2 
Si. oY eae 
a e 6 : (523) 
Z 








meererore 


ul 

Bo Sy V3(X9-X1) 

e = ee : (505) 
2 


iars can be substituted into Eq. (5.7) to give 
= 
d. 











Towa = simies 
Z 
and , | ie 
7 Sy XQ7-*] 
“no i) = |S 
Z 
Since these are equivalent, the estimate for S_ is 
2 “Lo 
x ,~X x,-X 
7 aaa 2a 
Ss Se SS (5.10) 
The standard deviation for So is given by 
Z 2 2 2 as 
: Xo AS, x) AS, 2 ; 
So = Soi | ge) ee | ee “5, Ge eee 
Ze ail 1 a 2 


The transfer function of the source is obtained with a simple 


vis 10n, 


Saye] 
Z.(V;) = 0.) (Sale, 
whereby 5 > 
ASV ) Z ARN al) 
AZ (V,) = Say Q i At(v.)° . (5.03) 


Pee hECT OF THE DISCHARGE 
In order to make a valid comparison between experimental results 
and theoretical predictions, the predictions should be based on 


values which represent the experiment. In this case, to compare 
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the Ingard-Schulz predictions with experiment, the effects of 
longitudinal and radial thermal gradients must be taken into 
consideration. Since the electron temperature and electron density 
are not known, they must be estimated. Also it Should be noted 
that the discharge does not fill the region between the two diodes 
but, using the values from Table 1, occupies only 55.6% of the 
Space. When there is no discharge, this does not affect results 
but when the discharge is on, the region is inhomogenous and the 
attenuation constant obtained is the average value for the entire 
Space between the diodes. 

Several assumptions must be made in order to predict the 
experimental results. The first is that there is no transition 
Peperane thewerds cf the disichangersregien. Thais Mmeams thet the 
predicted value for the effective attenuation constant should be 
the average of the value for no discharge and that in the discharge, 
where the values are weighted by the lengths of the two regions. 
Due to the nature of the interface between the regions, reflections 
are considered negligible. Temperature measurements taken between 
the discharge and diode 2 support this assumption. They showed only 
a one-degree higher temperature than measurements far removed from 
the discharge. The two diodes are separated by a distance d, and 


O 


the main discharge elements by d Sie eae denote 


hs on, ES Gos 


the attenuation constant in the discharge and outside of the dis- 


charge respectively, the effective attenuations constant Po 1S given 
d d,, 


dk 
eowh  SL 
O dy, Ol diy OZ 


by 
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Other assumptions concern the effects of the discharge through 
thermal gradients. During their investigation of the neutral 
temperature in medium pressure discharges, Hayess, et.al. [49] 
observed a radial temperature gradient under all discharge con- 
ditions. This will affect the propagating sound. Hayess [50] 
considered the problem in conjunction with his experiment. He 
concluded that for frequencies below the first radial resonant 
frequency, only the fundamental mode, a plane wave would pro- 
pagate. He further concluded that the temperature difference 
would not affect the attenuation greatly since only its square 
root appears in the attenuation factor. The cut off frequency 
for this experiment is about 9 kHz. The discharge region can 
Pemenoracterized by two temperaéurt@s for two Gititezent purecces. 
The longitudinal gradient in the tube affects the pressure in 
the tube. For this purpose, the dischaige can be characterized 


by its mean temperature T Since the temperature outside of 


M~ 


the discharge, ee did not change, the tube pressure with the 


discharge on is given by 


oe oe ee oe. 
ee wa Po qd, i To 
] - ae 1 - 7 1 - aati - = | 
1 M / M 


where d,. 1s the length of the tube. The uncertainty in the 


pressure 1S given by 


OT 0 


2 2 
t t 
ap? = op! AT 2 + opt Ap 2 : 
M M Po O 
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(seb 
’ 
Since <1 = 4054 = , this can be expressed as 
O 
M Ty 
V2 Z Z 
( se"\". ly on. og 0) 
f e 
7 1-2(.-z)™ Ml no 
do Tu 


Under the assumption that boundary effects are the dominant 


loss mechanism, the temperature at the wall, T TSeiMeOr tants Lor 


Ww? 
the discharge region. This 1S because these losses Bceue within 
two or three mean free paths of the boundary. 

To calculate the predicted attenuation constant with the 
Ingar d-Schulz equations, the experimental parameters must be 
converted to the variables used in the computer program: %, Noes 
Nee Te? and T° The wali teiperature can be used for tn: The 
neutral density can be calculated from p' with Ty? using the 
expression 


es CRO OU sxc 107° pee 


eae 
where p is in torr and NS in cm ~. The relative uncertainty in 


the density is given by 
AN \2 \Z2 49 2 
- P 
Since they were not meaSured, estimates must be used for the 
electron temperature and the electron density. With the similarity 


parameter Rp as an entry, ‘ips may be obtained from von Engle [40] 


with an estimated uncertainty of 10%. Estimating the electron 
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density requires a more complicated process. Venzke [51] has 
plots which give E/p as a function of i/R with Rp as a curve 
parameter. The value of E/p thus obtained can be used to obtain 
meee electron drift velocity from page 525 of Mdbaniel eons With 
the assumption of a uniform current density across the tube, the 
electron density is given by 


I 
ee, 
PARLE 
: 6 Dde 


Of course it 1S known that the current density 1s not uniform 


across the tube but neither are other factors Such as neutral 


density. The entire acoustic problem is treated in a one-dimensional 
sense. Because of the approximations involved, the electron den- 
Beas AaSSumed to be uncertain by 502. 


mae relation of Boe the attenuation constant, to 8, the 
attenuation co-efficient, can be viewed in two different ways. 
rea. = Bo Vw , the two ways are equivalent. One way is to write 
= = ll). 
B= p(w = 1) 
The other way is to write 


B= B(W)/ Ye. 


The latter form is preferred since a departure from the assumption 
of boundary loss domination will be readily apparent from a plot 
of B as a function of WwW. It should be noted that Figure 2, which 


B.-P 
ieee plot of sree , would not be changed by switching B Ole 


Pur 


Since Pa should have no dependence on w, the value of W should not 


be important. To be certain, the value corresponding to the nighest 
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and lowest experimental values of W should be used and the 
dependence checked. 

Since al]. of the input quantities are estimates, with an 
uncertainty which can be treated as a standard deviation, the 
estimate of A obtained will be uncertain. The uncertainty in By 


is given by 





2 2 D , 2 2 
eo . OT n OT ‘ e ON. : e ON. ss n ou) ‘ 


The partial derivatives can be obtained with the computer program 


by the expression 


= -|- = —_ . 
Ee = one os G8) ae S) 
Ox. AX. ; 


If AB 1s particularly sensitive to any of the variables, the 
process of evaluation should be reviewed to see if the uncertainty 


Gare be LYeduced. 
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VI. TREATMENT OF DATA 


The choice of data handling procedure was strongly influenced 
by problems associated with the experiment. The strong reflection 
from the tube ends was one of the most important factors in all 
considerations. The need for gated signals rather than continuous 
waves , in order to study the effect of the discharge, meant that 
the data could not be taken with meters. The multi-component 
spectrum associated with the bias for the gated sine waves meant 
that the signal had to be Fourier analyzed. 

Due to the strong reflections, the signals could not be 
analyzed in the continuous mode with a Correlation Function 
Computer such as the PAR 1OOA. When the gated mode was used, the 
time base was limited to the period from the triggering signal to 
the start of the reflected signal. Since this time base was 
sampled only 100 times, the resolution was very poor. Preliminary 
analysis of the signals indicated that the source transfer function 
decreased at low frequencies. Since this could help to check 
theories of the production mechanism, the low frequency components 
were of particular interest and resolution was needed. 

A Boxcar integrator was also used to recover the signal ona 
scale compatible with transcription inte numbers. However like 
the Correlator, this utilized only one signal waveform at a time 
and.it wasted most of that. In addition, the data had to be con- 


verted into numbers, and card punched for computer processing. 
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Even using digital voltmeters to obtain the data, obtaining more 
than one hundred twelve points from a waveform was very difficult. 

For these reasons, the process of recording on tape and 
direct conversion to digital form was finally selected. [It offered 
the greatest resolution, processing possibilities and error re- 


Bewiery capability. 


A. DATA REDUCTION PROCEDURE 
1. Analog to Digital Conversion 

After completing the data recording, the data tapes were 
reviewed with an oscilloscope. One minute segments were cut from 
each of the original tapes and spliced together. They were 
separated by long sections of white leader tape for easy identifi- 
cation of the segments. 

A hybrid computer composed of a COMCOR 5000 analog computer 
interfaced with an SDS 9300 digital computer was used for the con- 
version of the analog data to digital form. The analog side of the 
hybrid computer was set up to receive the analog signal and provide 
the appropriate signals to the digital portion for digitizing the 
information. The gain of the data channels was 50. A channel 
corresponds to the waveform recorded on it. Logic boards were 
used to start sampling at a frequency of 10 kHz when the trig- 
ger Signal reached a variable trigger level. Two switches were 
used to start and stop the conversion process. The number of 
samples taken following each trigger Signal was controlled by the 


digital program. After converting 1200 samples of each of the 
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three data channels to digital form, the information was written 
on a magnetic tape. The program used for the conversion is a 
standard program available to all users of the hybrid computer. 
The parameters for the program are the number of channels to be 
sampled and the number of samples per channel to be made for each 
record. The logic involved in the sampling process is different 
for most conversions and this required a slight amount of patching 
on the logic patch boards. 

Ze oeven to Nine Track Conversions 

The information written on the 7 track magnetic tape by 

the hybrid en is in the form of 24 bit binary words. To be 
processed by the IBM 360, these must be converted to 32 bit binary 
Wemge. the only programming invoived in (his is @ prcogwan te aead 
mnemystrack tape, call a SeuaEmeeumccite and write the information 
on another tape, usually 9 track for increased efficiency, for iater 
processing. For this job, this step was also used to separate the 
samples by channel. 

&. Ensemble Averaging wath Data Rejection 

The diode signals contained noise. The noise was random and 

consisted chiefly of positive-going spikes. A low pass filter could 
not eliminate them without distorting a significant wOrtion Of Ehe 
acoustic waveform. In addition there were a few bad trigger signals 
recorded. Thus it was known from the emery. that some of the 
data records were bad and should not be included in the processing 
but it was not known which they were. Since the entire data set 


consisted of ten files, each of which contained from 80 to 
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records of 3600 words, screening became a major task. The 
following procedure was finally found to be most efficient and 
effective. 

Each file corresponded to an experimental run and all 
records in that file represented the same waveforms plus noise. 
The records jin a file were all averaged together and the standard 
deviation for each channel was obtained. Following this, all 
records were compared with the average. Those with less than 75% 
of their points within 2.8 standard deviations of the average were 
dropped from further processing. Using only the acceptable records, 
a new ensemble average and standard deviation was obtained for 
each channel. These were used in a final ensemble averaging pro- 
cess in which only those points witnin 2 standard deviaiicis oF 
the average were used. This provided an ensemble average of the 
waveform and the standard deviation of the points involved. 

4. Delay, Overlay and Zero 

At this point, the data consisted of 1200 samples for each 
of the three waveforms, representing their average values, and the 
standard deviation of the waveforms. For the sampling process 
1200 samples were taken so that all of the desired waveform of the 
Signal at point 2 was sampled. This meant an undesired signal was 
included in the samples of the signal at 1, namely the reflected 
Signal. In order to make a valid comparison of spectral components, 
the time samples of the waveform which is transformed must contain 


the same number of samples. Reducing the number of samples of s, to 
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those necessary to avoid the wflections in s, means the starting 


1 
point for the transformation must be delayed. Since phase infor- 
mation was not going to be considered in the analysis, and time 
delay merely changes the phase of the components, delaying the 
Signals would not affect the results. Therefore the signals Sy 
and S, were delayed by choosing a starting point for each such 
that 600 sampies included the waveform up to the point where the 
reflection started. 

The diodes interacted electrostatically with discharge. 
This caused a small coherent square wave in the diode output 
which coincided with the time of the current pulse. fees inclusion 
in the waveform to be transformed would give erroneous spectral 
.1 ere Se 


components for the acoustic sign Delaying weee cuffeac 


eliminate this unwanted signal from s in order sto elimina ce at 


2° 


from s without introducing unwanted spectral components, the 


1? 
values of the samples in the region containing this unwanted 
Signal were replaced with values taken at random from the adjacent 
region. Since this adjacent region contained only samples of the 
waveform before arrival of the acoustic signal, they only repre- 
sented background noise. | 

Since the initial portions of the signals S) and So 
represent no signal, i.e. the period before arrival of the acoustic 
Signal from the sound source, itS average value should be zero. 
mere ionais were offset to zero by first determining the average 


value of the no signal region and then subtracting it from ali 
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values. For the current waveform, this region followed the 
current pulse. The three operations of delaying the signals, 
overlaying the electrostatic signal and establishing the zero 


level were performed in one program. 


fae OPERATIONS ON THE TRANSFORM 

Two sequential programs were used totransform the waverorms 
and obtain acoustic informations from them. The first program 
took the input waveforms, corrected them for gain and diode 
sensitivity and determined their spectra. The second used the 
spectra to determine the transfer functions of the discharge 
region and the source. 

Input data for the transform program were the three ensemble- 
averaged waveforms, i(nT), S,(nT), S5(nT) and their standard 
and the diode sensitivities. 


Geviations, Ai, AS AS 


Two uncalibrated amplifiers were used in amplifying and 


le a 


recording the diode signals. These were calibrated as part of the 
processing. While recording data, both amplifiers were fed a 
common Signal consisting of 2 millivolt peak to peak, 5000 Hz, 
tone burst. This calibration signal was processed in the same 
manner as the other signals. In the first part of the transform 
program, the spectrum of this signal, in the diode channels, was 
Fourler analyzed. Then the spectral component closest to 5000 Hz 
(5040) was used to calibrate the gain for each diode channels so 


that this component was equal to one. 
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The diode sensitivities, Supplied as inputs, were calculated 
earlier from experimental data and expressed in terms of milli- 
volts/torr. The sensitivity of the current detector was 
Dx 10°? volts/milliamp. The sensitivities and corrected gain 
factors were used to convert the values into the correct magni- 
tudes of the physical quantities involved, milliamps for the 
current and torr for the acoustic signals. 

The transforms of the corrected waveforms were then calcu- 
lated using a computer library subroutine "RHARM". This program 
requires that the number of points be an integer power of two. 
There were 600 points available and the nearest power of two was 
512, leaving 88 excess points. This means that there is some 
Wm@eelicalnty aS tO where to start the analysis. or *#hie reason, 
the waveforms were analyzed five times, each time advancing the 
starting position by 10 places. The average value and variance 
of each component was calculated from the five spectra. 

The output of the transform program was the average spectrum 
and the variance of its components. Only the first ‘twenty five 
of these were passed to the next program. The magnitudes of the 
remaining components were less than 5% of the fundamental and were 
characterized by a great deal of scatter. 

The next program used this data, and calculated the acoustic and 
source transfer functions. The equations of Chapter V were used 
for the calculations. A computer library subroutine "LSQPL2" was 


used to obtain the least squares determination of Bo: For this 
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calculation, the data was weighted inversely proportionally to 
the Square of its uncertainty, with weights normalized such that 
their sum was equal to the number of points. 

For the data corresponding to no discharge, Eq. (5.10) was 
used to obtain an estimate for the sound spectra at the source. 
This was divided by the sound source current spectra to obtain 


fimestransfer function of the sound source. 
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VII. PRESENTATION AND DISCUSSION OF RESULTS 


A. AVERAGE WAVEFORMS 

The averaging and screening method worked very well although 
the program could be rewritten to perform the job more efficiently. 
Miyeethree waveforms for the discharge-off case at 7.65 torr are 
shown as Figures 10, 11 and 12 just before being transformed. For 
contrast, the first 2/3 of the average waveform for Sy for the 
same run is shown as Figure 13 with the same scale. The second 
puise is the reflected signal from the end close toa the sound 
source. This shows the effect of the delaying, overlay and zeroing 
operations. Unless otherwise specified, the remainder of the 
es cs displayed in this chapter were derived from these wave- 


forms. 


BemeicOUSTIC TRANSFER FUNCTION 

The points of the acoustic transfer function showed a rea- 
sonable amount of scatter about the predicted form. The greatest 
degree of scatter was found to be associated with those components 
whose contributions to the signal were small. Figure 14 shows 
the first twenty five spectral components of the sound at point 2. 
Figure 15 shows the transfer function ie a linear form. Figure 16 
is the same information expressed as attenuation of sound in 
decibels per meter as a function of frequency. The curve repre- 
sents the least squares fit to the data under the assumption that 
boundary losses dominate. The goodness of fit found in all cases 


eupeortsS the validity or this assumption. 
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Figure 10. Sound Source Dischaxge Current. P = Jaos Orr . 
Main Discharge off. Time Scale: 2.048 msec/in. 


Ver t@ical scale 1s axbitrary. 


Froere 1%" Sound Signal at Peint Ieeidelayed, overlayed 

and Zer0ed);. P =87.,058ceru9™ Naan discharce 
off. Time S€ele. 2.0456 maeec/in. Vertical 
scale 15 arbitrary. 
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Figure 12: Sound Signal at Point 22, (delayvea, overtaved 
and zeroed). P = 7.65 torr. Main discharge 
Off. | Dime Scale 25040 msice7 ine vcr tlc an 
SCalec iS arbitrary. 


Re pe 


prgutre 14. Sound Signallar Powe Ie (unecirectaed) ; 
P = 7.65 torr. Main discharge off. Time 
Scale = 2.046 meee/In. Vertical scale is 
the same as Figure ll. 
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Spectral Decomposition of Sound Signal at Fo 
This 1s the-sayeeoersl iver i905 itton © 
form SHOWN Tyan eres eas eg lane 
simple connection of points. 
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Pagitiecr 1:5: S,(V)/S (Vv) in Dependence on v*? for a Pressure 


of 7.65 torr, “Main discharge Off. Breer bars 
represent One standard deviation. Curve 
represents the weighted least squares fit to 
the data. 
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Attenuation in Dependence on Frequency for a 
Pressure of 7,65 torr. Main discharge off. 
Curve represents the least squares fit to the 
Gite, This fare ormesponds to figure 15. 
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C. SOURCE TRANSFER FUNCTION 

The source fyvapeeennees was found to have the low fre- 
quency drop off indicated by the preliminary observations. At 
frequencies above a pressure-dependent transition frequency, the 
frequency characteristic did follow the predicted form. The 
magnitude was also found to agree with the predicted values. The 
frenster functions are displayed as Figures 17, 16 and 19 for the 
femeeing pressures§7.65 torr, 12.7 torr and 16 torr respectively. 
The structure is believed to be due to scatter in the data since 
the points with the greatest departure from the smooth curve 
correspond to those spectral components with the largest relative 
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Two aspects of the results can be used to test the theory 
of sound production from a discharge. One is the frequency de- 
pendence and the other is the magnitude. 

a. Frequency Dependence 

The low-frequency drop off is contrary to the fre- 
quency dependence predicted by Schulz's results. This effect is 
believed to be due to neutral particle diffusion under the in- 
fluence of thermal gradients. It is a low-frequency , low-pressure 
effect which should appear only when thermal energy sources provide 
the acoustic energy. This will be Are cueeee later. 

Recording to Schulz 's predivetiton tie etranster function 
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Figure 17. 
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Source Transfer Function in Dependence on 
Frequency for p = 7.65 torr. Main discharge 
Off. Vertical seale ai ints of (newt on/n< ) 
per ampere. Solid curve represents predicted 
form fitted at one point... 
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Figure 18. 
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Source Transfer Function in Dependence on 
Frequency fer p = V277 stouces Main discharge 
Off. Ver ticaleces te ities. Of (newtons/m~) 
per ampere. Solid curve represents predicted 
form fitted at ene poi ce 
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2 3 so 5 
YJ (kHz) 
Source Transfer Function in Dependence on 
Frequency for p = 16 torr. Main discharge off. 
Vertical scale in units of (newton/m“) per 
ampere. Solid curve represents predicted form 


fitted at one woman. 
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This prediction was tested by fitting a curve given by 





Cc 
_ 1 ae 
Y > eee oe) (Joely) 
through the transfer function data points. The coefficient C 


it 


was evaluated by causing the curve to pass through the point at 
vy = 2450 Hz. 

Schulz's derivation was based on planar electrodes 
separated a distance "a". The experiment used cylindrical elec- 
trodes of length 1% cm separated by a distance of 3 cm. Therefore 
the effective length of the discharge could lie anywhere between 
3 and 6 cm. The value of "a" was adjusted within these limits to 
@@tain a curve which fitted the data reasonably well. For 7.65 
torr, the value of a was 3.5 cm: for 16 torr, the value was 
4.6 cm. Since the pressure affects the discharge, there was a 
Pessibllity that the value of "a" for 12.7 torr could be inter- 
polated from the other two values. The value obtained by inter- 
mobration was 3.93. This gave a good fit to the data. The smooth 


curves shown in the figures are the fits of equation (7.1) to the 


data. 
b. Magnitude 
Using Schulz's equation, the source transfer function 
1s predicted as sin(=— 
Z. = (y-1)H 5S 
= Iw 
The energy transfer rate from the electrons H is the 
Same as used earlier in conjunction with wave amplification. [It 
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may be expressed in terms of the characteristic time for electron 


heating as 2 


H = a 
Tye) 


Zz : 3 

The calculated value of c /T was obtained with the 
computer program used for evaluating the Ingard-Schulz equations. 
This was accomplished for the transfer function at vy = 2540 Hz 
and p = 16 torr. The discharge current spectral component at 
this frequency was 6.4 ma which was estimated to correspond to 

: 10 a ae 

an electron density of 2 x 10 cm ~. In determining this density, 
the cross sectional area of the electrodes was used since the dis- 


charge did not extend, in the radial direction, significantly be- 


yond them. The procedure used was the same as outlined in Sec. V.E. 


The predicted value of Z Senne wa S 


S 5 
Ace ampere 


The corresponding value, experimentally obtained, was 


Z 
7 : nt /m 
cone ee ampere 


The agreement of the two values is good. Considering 
the number of approximations made, order of magnitude agreement 
would have been considered satisfactory. No error analysis was 
conducted on this computation for that reason. 

Except for the departure at low frequencies, the re- 
sults are consistent with the theoretical predictions based on the 


assumption of energy transfer by the electrons in the positive 
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column. There are aspects of the problem which were not examined 
by this experiment but which should be studied to gain complete 
understanding and to understand other esperimental results. A 
more complete investigation would include phase as well as magni- 
imgde information in the transfer function. In addition, the 
source of the sound could be localized and momentum transfer 
aspects investigated. Although this experiment indicates that 
electron heating of the neutrals in the positive column provides 
a mechanism for sound production in glow discharges, it has not 
discounted the roll of ions in sound production under other 


circumstances. 


im Low frequency Departure 


During preilminary wOrK in COnjumeti0n With this eeperi- 
ment, a large number of acoustic waveforms were observed. One 


characteristic of the sound signal produced by a pulse was noticed, 
which is reflected in the drop-off of the transfer function at 
low frequencies. Figure 20 shows part of the waveform received at 
point 1 during one of the runs. Me to attenuation of the high 
frequency components, the features are not as sharp as some which 
were observed at closer positions during earlier work. The main 
features are visible however. 

The first and most significant feature is that the vertical 
distance from 1 to 2 is approximately equal to the distance from 
3 to 4, with the latter always slightly greater. The next feature 
is that the slope and curvature of the curve from 2 to 3 is approxi- 


} 3 


mately the negative of that for the curve from 4 to 5. Ine Gro 
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Figure 20. A Typical Sound Pulse Produced by a 
Rectangular Current Pulse. 
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curves suggest an exponential process with a long time constant. 
A last point is that point 5, where the slope is approximately 
zero, is lower than point 1, the level before the signal. The 
polarity of the figure is such that point 1 corresponds to the 
ambient pressure with upward deflections indicating an increase 
in pressure in the discharge region. 

The following description of events occurring in the dis- 
charge region during the creation of the sound pulse is postulated. 
When the discharge current increases, the energy transferred to 
the neutral gas increases, thereby increasing the temperature of 
the neutral particles. This temperature increase 1s si rapid that 
the netural particle density remains essentially constant. Asa 
C@nsequence, the pressuxrc changes almost directly with the tem- 
perature. The neutral particles move away from the discharge, 
radially as well as axially, under the temperature and pressure 
gradients. As the concentration decreases, the pressure decreases 
toward the ambient value. The pressure change, AP > which would 


take place because of the diffusion is 


4p, = KANT, 


where AN is the change in density. 

When the heating processes are removed, the neutral gas 
cools rapidly to the ambient temperature. Since the cooling 
process is much more rapid than the particle diffusion, the pres- 


sure at ambient temperature would be the equilibrium pressure 


Nie Ws 





minus the pressure difference caused by the density change. This 


pressure difference, AP5» 1s given by 


AP, = KANT. 


>< 


Since Ty To 5 APy ~ Ay The neutral particles would then 
diffuse back in, to return the region to the ambient density. 

The diffusion process is slow, with a long time constant. 
Under the action of a high frequency signal, there would not be 
time for a significant amount to occur between cycles. Coneeenen tye 
the acoustic pressure would not reflect any density changes. At 
the other extreme, for very low frequency temperature variations, 
the pressure eee would be negligible but the density changes 
would be large. Therefore it can be predicted that the effect 
will become appreciabie at low frequencies. 

Notice that in the os uae process, no assumptions were 
made concerning the heating process except that it is spatially 
uniform in one region and that particles are free to move in and 
out of the region. Thus the assumption of one dimensionality for 
such a process occurring in a tube could be a valid assumption. 
Therefore Ganniele assumptions will be considered valid and the 
failure of the prediction at low frequency attributed to the 
wave equation. The wave equation Modis (oes it does not take 
diffusion into consideration. 

TAyditrusion is the important faccer sin the discrepancy ar 
ioweeeequencies, then the effects will become amportant for fre- 
quencies whose period is of the order of or longer than a character- 


fmtice diffusicn time x This time 1s the time for partic.e 
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diffusion for a mean displacement of some characteristic dimension 
of the discharge region. Therefore, the characteristic frequency 


1S given by McDaniel [39] as 


it D 
Vv = —_—_ => —_— 
Cc is Z 
C 
where € is the characteristic dimension. The effective diffusion 


coefficient D is given by Present [36] as D -% where v is the 
mean velocity and £ is the mean free path. 

Since the mean free path varies inversely with pressure, 
this predicts that the frequency at which the effect occurs will 
vary inversely wach pressure. An examination of Figures 15, 16 
and 17 shows that this is the case: at 16 torr, vee SOO Wi saat 


Mort 5 WEeO 
Cc 


tr TT ~ +h _ 
v Lia 3 vile r~ Sac t 


WN 


ems torr, My = 750 antiga. Cc 
Dee being roughly constant as it wou for the predicted dependence 
On density. However, substitution of the discharge length "a" for 
the characteristic dimension gives a frequency which is three ordeis 
of magnitude too low. Thus, although diffusion effects may be 
responsible for the low frequency effects, the simple argument 
used to determine the critical frequency did not predict the 
Cemreect result. 

This diffusion effect 15 a new problem for acoustics. ist 
is a low frequency, low pressure effect and will probably only 
Seeur in conjuction with a thermal acoustic source. The excitation 
of sound by a heat or energy source is self-rectifying since the 


energy transfer into the system cannot be balanced by the removal 
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of heat or energy from the system by the same process. This 
means that the energy source term H in equation (2.1) can not be 
represented by the form 


H( ryt nee 


Since such a representation implies that energy can be removed by 
the source with the same efficiency that it transfers energy into 
the source. 
| To predict the correct low frequency characteristics of 

sound produced by a discharge plasma, the derivation of the wave 
equation must take diffusion and heating effects into consideration. 
The results of the derivation should be a differential equation 
capable of predicting the neutral density and temperature when 
See ca fcr the zero frecuency case. Such a derivation 1s con- 


Sidered beyond the scope of this investigation. 


Peer PECT OF THE DISCHARGE 
1. Predicted Attenuation 
When the measured values of axial temperature were com- 
pared with the predictions of Venzke et.al. Come the agreement 
Wesewithin /%. This verification of the predictions of that paper 
was considered sufficient justification to use its predictions for 


T._, the wall temperature and T 


Ww M? the mean temperature with the 


same relative error. These values were substituted into the ex- 
Pressions.of Section V.E. to obtain the neutral density in the two 
regions for each run. The steps outlined in that section were 


motlowea to obtain estimates for we and T . With the density. 
e 
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wall temperature, electron temperature and electron density as 
inputs, the attenuation constant was calculated using the Ingard- 
Schulz equations for W = 31,400 and w = 600. 

The results obtained were that the zero ionization value 
of the attenuation constant was essentially the same as the value 
predicted for a discharge under the experimental conditions. The 
value of B., calculated in the manner mentioned earlier was larger 
at the higher frequency, which indicated that bulk losses were not 
completely negligible. The largest difference between the values 
one B for the two frequencies occurred at 7.65 torr where the 
difference was slightly less than 3% for the ae ioneed’ajanneee The 


greatest lonization effect predicted was a .5% reduction in 


feeeennatcion for the stwo dischare@m runs made with @ filling pmes- 
Sire Of 7.65 torr. In other words, the Ingard-Schulz equations 


predicted that no plasma effects would be discernable under the 
experimental conditions. 
For the plasma associated with a filling pressure of 


7.65 torr, the electron temperature is was estimated to be 


l ll +3 


fees 10 K and the electron density was 5 x 10 cm ~. With these 
: | _ | Leese 

values and the values for the neutrals, NES Zee 1O@oeem sand 
Th = 323K, the values of some of the plasma parameters were cal- 
culated. The characteristic collision frequencies are; "Wo 29 x107 
and [US pon 10/ sec. The electron and ion plasma 

; 10 §  -1 
frequencies (W,, W,), are 400 = TO] ance S = f0eucn 5 The mean 
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free path of a neutral particle is 7.0 x 10 *° cm. While these 
values are different for the other discharges, the values are 
representative of the orders of magnitude of the various para- 
meters. 
2. Comparison of Predicted and Experimental Values 

The assumption of boundary loss domination was made 
earlier. Since the failure of this assumption would occur at 
higher frequencies, where the data were weighted the least due 
tO their larger relativeefror, the fit to the data should reflect 
the low frequency characteristics. This was tested at the value 
Or BY predicted for W = 600. The predicted values are shown in 
Table [II under the column marked Boe The experimentally deter- 
eeee volwe cf BG and its standard deviation are shown in the 
columns marked Bo and AB O° 

The ratio of the difference, between the predicted and 
experimental value of Bo? to ene Fauceed deviation of Bo 1s 
shown in the column headed . The magnitude of these 

O 

numbers are all less than unity, indicating that the differences 
between the predicted and the experimental value are all within 
one standard deviation of the result. The greatest difference 
between the predicted and experimental values is 13%. 

However the signs of the difference are correlated. The 
second column, headed "polarity", indicates the polarity of the 


discharge with respect to sound propagation. A plus sign means 


sound propagation in the direction of electron drift motion. A 
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zero is used to indicate the discharge-off state and a minus sign 
indicates sound propagation counter to electron drift. The sign 
of the differences is positive for the no-discharge cases. For 
the discharge cases, it is negative. The negative sign indicates 
that the predicted attenuation is greater than the observed 
attenuation. 

SB, (Discussionvof the Results of theyGonparisoen 

The results show that for the 1o0nized cases, the predicted 
value of the attenuation constant is greater than the observed 
value. Although the two values lie within one standard deviation 
of the experimental value, the consistency of the reste indicates 
that it is due to either a systematic error or to a plasma effect. 
Miemeeompernture dependence ef the attenuation caonstamt at wanstant 
pressure is such that if the wall temperature used for the com- 
putation were too low, the predicted values of B would be too 
high. The three parameters which influence the predicted values 
to any degree are W, Ne and To: The error must be due to either 
Ne or cha Since although the value of Wis the same for the ionized 
and un-10nized cases, the results of the comparison show that the 
sense of the errors involved is different. 

In calculating the predicted values, the propagation of 
error was followed. Although the mean temperature in the dis- 
charge was assumed to be uncertain by 10%, it did not affect the 
error of the corrected pressure significantly. That error was 


affected primarily by the uncertainty in the value of the filling 


124 


pressure, which was estimated to be 5%. If the error in filling 
pressure was Systematic, it would have affected all cases in the 
same manner. Therefore neither an error in the mean temperature 
nor an error in the filling pressure is believed to be a source 
of systematic error. 

The only remaining source of systematic error would be the 
value of wal] temperature used. This affects the results through 
the density, calculated using T. and pp’) andmdarectly ..elhe cuaor 
in the wall temperature was conservatively over-estimated to be 
ne , which means a 25% error with respect to the increase above 


ambient temperature but only a 3% error with respect to the ab- 


solute temperature. Representing the error in ten due to the 
ermor i: - by Op. ATOM eS aM © pla SO LORY 
oe) 
T ON | or OT (1) 
Simce 
N = aus 
i ict 


Bor constant pressure, 


which gives 


Using this expression, 5B was calculated for es LO. ~ Lhe 
results are displayed in Table III. It is apparent that these 


values are not large enough. to account for the discrepancy between 


results except for runs 7 and 9. It should be noted that if the 
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value of the filling pressure used for predicting Bo fotet bins. 7. 
8 and 9, were high, then the correct value for the prediction would 
be higher for all three cases. This would decrease the dis- 
crepancy for the no discharge case and cause the results for the 
ionized cases to be more like the others. The conclusions to 
pemdrawn from this discussion are that while a Systematic exror in 
the value of the wall temperature used for predicting the results 
1s possible, it is not likely to be the cause of the difference 
between the experimental and predicted results. 

The two runs where this error should be the smallest are 
runs 2 and 3, Since the temperature increase is less at lower 
pressures. The difference for these two runs is not significantly 


~aenrsa @ 


imemer. if the results for Liliing pressures 7.65 and 16 torr Sine 


CN 


are compared, then it appears that the polarity does have an effeci. 
However, there is insufficient replication to support such a con- 
Glusion. 

The conclusion supported by the results is that for these 
experimental conditions, amplification of the sound may have 
occurred and a polarity effect may be present. The observed effects 


are within the experimental error. 


Pee VALUAIION OF THE YATSUI EXPERIMENT 

WatSul, et.al. [8,9] conducted their experiment in Argon at 
low pressures and high frequencies. For a comparison of their 
experimental results with the Ingard-Schulz prediction, the 


Ingard-Schulz equations have now been evaluated with parameter 
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values corresponding to two points in the YatSui experiment. For 
0°? Hz and I = 5 mA, they reported maximum amplification 
Sceurred at 1.5 x Tome torr and damping occurred for pressures 
greater than .31 torr. For the condition of maximum amplification, 
the ratio of the received signals with and without a discharge was 
3.3. The predicted value for this set of experimental parameters 
will be discussed first. 


Since the Signal traveled only a fraction € of its total 


travel in the discharge, the ratio of the two Signals is 


Pot (ier =e 2 
' (CB \+(1-e)p)d — -£ (8 -B) a 
CSie: = & See ee ome, e 
Ve a7 bd 
where dis the total distance between source and receiver. There- 


fore, the difference, Pye between the two attenuation coefficients 


can be determined from their data by 


iemaenominator, €d, 1s the distance traveled in the discharge, 
38 cm. Therefore the experimentally determined value is 


B28 = Se in 343 =e 2 nore aa 


Using the procedures outlined in Section V.E., the electron 


& 


temperature was estimated to be 5.5 x 10 K and the electron 


10 -3 


density was estimated to be 5.1 x 10 cm ~. With these values, 


6 - 
and WwW = 6.28 x 10°, No = 5x 107 cm 3 and TY = 300K; the quantity 
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Ee was computed using the computer program. This calculation 
gave Poe oer er, which 1s lower than the exper imental 
value by two orders of magnitude. USing liberal estimates of the 
errors in the input variables, the uncertainty of this result is 
Oe, ie Sys 107°, Thise1sS not enough to ive ordermcr 
magnitude agreement. 


The Ingard-Schulz predictions were also evaluated for the 


discharge conditions corresponding to the critical pressure, 


: - - 0 
Sol ctOrr, uSing: Ne = 1.1 x 1077 cm 3 No = ex 107 ; 
T. mt .O X 10" and ue = 300K. The value predicted was found to 
-10 2 
De B -f = - 9.54 x 10 : cm I ; 


p 


Since this difference would hardly be detectable, it is in 
agreement with the yYatsui results. Further agreement with the 
results is the prediction that,for pressures above this, damping 
Occurs. 

Thus, although the Ingard-Schulz equations do not predict a 
G@irectional effect, its predictions are, to a certain degree, in 


agreement with experiment. 
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VIII. SUMMARY AND CONCLUSIONS 


This study has examined the creation and modification of 
sound by a discharge plasma both theoretically and experimentally. 
Following a review of earlier experiments on the creation of sound, 
a mechanism for sound production by ions was proposed and shown 
to be feasible. Then a numerical analysis of the equations of 
Ingard and Schulz was performed. This analysis showed that net 
amplification would not occur. It also showed that the approxi-~ 
mate form of the dispersion relation was not useful because the 
criterion for its applicability was invalid. 

Incidental to the experimental study, the theory of the 
mobility-limited thermionic diode used as a microphone was de- 
veloped beyond its previous, elementary form. This analysis and 
the subsequent successful usage as an acoustically transparent 
microphone for two point measurements has increased its usefulness 
as an experimental tool. 

The diodes were used in an experimental investigation of the 
modification and creation of sound at medium ec Four ier 
transforms and computer data processing techniques were used in 
the analysis. of the data. The investigation into the modification 
effect of the discharge was ina different region of current, 
pressure and frequency from earlier experiments. The results in- 
dicate that an amplification mechanism may be operative and that 
there ce directional seftects Since tne Ingar d-Schulz 


equations, evaluated for the experimental conditions, predicted 
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no effect would be detectable, these results do not support the 
equations. The Ingard-Schulz equations were also evaluated for 
the conditions of the Yatsui experiment. The predicted values 
were not found to be in agreement with the experimental results 
but the trend with pressure was in agreement. The conclusion 
drawn was that the theory needed more development. 

The major contribution of this study is believed to be the 
information gained from what was originally the secondary experi- 
ment, the conversion of discharge current fluctuations into sound. 
This is the first known investigation in which the current magni- 
tude was SoueAee to the sound it produced. The results were in 
reasonable agreement with the Ingard formulation of the problem. 

A low frequency effect was observed which was interpreted as the 
mist Observation of the effects of particle diffusion in acoustics. 
The general agreement of Ingards predictions in both areas 

indicates that electron heating of the neutrals is the dominant 
source of acoustic energy. However several experiments have 

raised questions not answered with the present theory. One of 
these is the directional perce on discharge modification. Two 
others, associated with sound production, are the effects of high 
fields and the localization of the sound source near the cathode. 
These stand out as subjects of further experimental and theoretical 
development. The interpretation of the low frequency behavior of 
miemsound Source as a diffusion effect also warrants fur ther work, 


both experimental and theoretical. In addition, the study of the 
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Ingard-Schulz equations indicates that the increased attenuation 
to be expected under certain conditions could be a fruitful area 


for investigation. 
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APPENDIX A 


THE MOBILITY-LIMITED THERMIONIC DIODE AS A MICROPHONE 


Standard microphones have limitations which hinder their 


usage when studying sound at very low pressures. This is parti- 
cularly true in conjunction with plasma-acoustic work. There are 
several reasons for this. A very important one is contamination 


from gassy components. Another consideration is that the micro- 
phones were designed for operation at atmospheric pressure. 
Operation at reduced pressures may require circuit modification 
Since, for condenser microphones, breakdown can occur. Chalupnik 
ema, | 52) reported changes in sensitivity and frequency response 
with pressure in two commonly used condenser microphones (Western 
Electric 640AA, Bruel and Kjaer type 4111). The lack of a sensi- 
tive sound detector without the problems outlined above has 
slowed plasma-acoustic studies. Recently Dayton et.al. [53] re- 
ported the use of a mobility-limited diode as a microphone without 
the drawbacks of conventional microphones at low pressures. How- 
ever heretofore the theory has not been developed beyond its most 
elementary form. 
1. Theory 

If the cathode is sufficently hot, so that the operation 
is not emission limited, the current flow in a thermionic diode, 
Operating in a high vacuum, is limited by the effects of space 
charge. The current density is given by the Langmuir-Childs law 


which says that the density is proportional to the cathode-plate 
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potential to the three-halves power. The constant of pro- 
portionality is dependent on the geomentry of the current-flow 
situation. The three halves power dependence assumes an equi- 
potential filament; a condition which is approximately true pro- 
vided the voltage drop across the filament is small compared to 
the cathode-anode potential. If this condition is not met, then 
the current density varies with the plate potential to a higher 
power, up to 5/2 [54], the actual value of which is determined by 
the definition of the cathode to anode potential. 
a. Mobility-limited Operation 

ae presence of a background gas of Beutel parti- 
cles, the electrons will suffer collisions and their average 
Memieert, Will therefore decrease. This affects the snace charge 
field and the current density. In the limiting case, where the 
anode-cathode separation is much larger than a mean free path, the 
velocity of the electrons is limited by their mobility in the gas. 
Cobine [551 derives an expression for the current density ina 
parallel-plate diode with a mobility-limited current: 


Ein fini 
3 


ies Oe Oe 1O ieapeen (ere is) 


where d is the electrode separation in cm and » is the mobility. 
The derivation assumed the electron velocity is proportional to 
the field, i.e. the mobility is field independent. This result 
means that the current 1s proportional to the square of the anode- 
cathode potential. A derivation for a cylindrical situation under 


the same assumptions yields the same dependence on plate voltage. 
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Although a general treatment has not been attempted, the similarity 
of the derivation to that of the Languir-Childs law suggests that 
the result is geometry independent. With the additional assumption 
that aye mobility WU is inversely proportional to pressure, the 


current ina mobility limited diode is given by 


2 
V 


1 (6% a e.s 


For weak fields, the mobility is not field-independent 
WO and the drift velocity is given by 


1 
2 


E 
(= ; : 
ves) (A.3) 
In such a case, it 1S easy to show that the diode current for 


parallel plate and cylindrical diodes is given by: 


Bye 
V 

1 

“5 


Pp 





ioe 


(A.4) 


The field-dependent, the field-independent and any 
intermediate case as well as the transition from space charge to 
mobility limited operation are included in the general form 


“kV Kk, 


f= op (5) 
where C, Ky and Ko are dependent on the situation and must be 
empirically determined. 

beans Sensitivity 
To determine the sensitivity of the diode as a micro- 


phone, consider the simple circuit shown in Figure 21. If the 


effect of the sound on the diode is observed by measuring the 








Figure «21. DileGder@irecuit. “A -simplescreeu1t Whicrmesnemc 
used to study sound with a mobility-limited 
thermionuve diode: 
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Figure 22, Diode Current in Dependence on Pressu'’« fcr 


Typical Diede in ArGon and Neon: 


sof, 


voltage across R, then the Sensitivity, S, Of he varedevcaw cua 


can be defined as: 


OW oe 
S=- [ee (A.6) 


Using equation (A.5), it is easy to show that this is given by: 


k V 
i DC 
p V 
De 
ile a V=¥ 
DC 


where Vig is the d.c. component of the output voltage. This 
expression shows that the sensitivity increases with the d.c. 
voltage across R, in an approximately linear relation, Unewe the 
output voltage becomes comparable with the supply voltage V. The 
Signal voltage V produced by an acoustic pressure wave Dp 1s there- 


fore given by 


ad 


We So. (A.8) 


c. Frequency Response 
Equation (A.7) is not correct for short wavelengths. 
Implicit in the discussion thus far has been the assumption that 
the finite size of the diode would not affect results. Actually 
the diode averages the pressure over its length. When fringing 
effects are ignored, the diode agen is proportional to its 
length. Therefore a diode of length 6% contributes a current 


61 given by 50 
61 =f i (A729) 


where I is the current for a diode of length £, considered uniform. 


The total current consists of the d.c. component plus the a.c. 





component caused by a pressure disturbance. The current fluct- 
Hatton 1S -oiven. by 


T==2p. (A.10) 


AIM 


1kx=-iWt 
e 


For a pressure wave given by p =D 


» Which is propagating 
Poethe axial direction of the diode, the total current £ 2s given 


by integrating equation (A.9) over the Length of the diode, wherein 


the real part of the pressure must be used. 


fs 
aes o) 
1 = } Z GOST sc ccs 
O 
ther efore: eS sin ki (A.11) 
R Po ee 


The frequency response of the didde is therefore es 
G. Temperaturve effect 

The equations pevelapad thus far indicate that the 
diodes can be eal eee in terms of absolute pressures by extra- 
melation from the zero frequency case, 1.e. the d.c. characteristic 
of the diode. This requires that the pressure be varied and the 
current measured as a function of plate voltage at the pressures 
used. The pressure measurement is usually made at some point 
removed from the diode. Since the filament is hot, the effects 
of the thermal gradient between the diode and the measuring point 
must be considered. 

Many properties of particles in gases which are fre- 
quently expressed as functions of pressure are in fact primarily 
dependent on the particle density. In such cases, temperature 


effects can be compensated for by correcting the pressure 


ia 








temperature. For example, mobility is given as a function of 


pressure even though it is a function of density. 
b = w(p) 


If the data is for pressures at temperature i and the temperature 


of the point of concern is T and the pressure is Po» then the 


elas 


effective pressure, Po? £Or “entry Inte thevdata 1s 


Now consider the effect of the temperature difference 


on the diode characteristics, starting from Gauss's law: 


yo =. (A.12) 


and the expression for current density 


I ONE 


Assume electrons moving ina gas under a weak field so that 


i 
ee 
y= Ga 


p 
where Cy LS @ cCOnStant peculiar to the oas. Consider cy Pindri cai 
symmetry and negligible end effects. Since the current I is a 


constant, and related to j by geometry, the density can be expressed 
fiaeterms of the current and other factors then Substituted into 


Eq. (A.12). The result can be written in cylindrical coordinates - 


as: 
ays 
eet | 2 ee ae 
ia de dr amx£ C, /1 av) % 
( sr ) 
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This can be re-expressed as: 


fe 
xr dvj* da fr av j_ he 
dr dr dr 211 


The pressure p can be expressed in terms of the 


(A.13) 


“5 [° wr 
H 


pressure, Pj; at the measuring point as 
ee tthe) ) gegtGs)) 


With this substitution, Eq. (A.16) can be integrated 


372 p.° ale 
( = -355 PL ery eter! 





1 
Therefore 
ee) 
2 
p i 2/3 
dr 2 em C, 16 
= 
F(x) -( fea yr ei | 
Oo 
mom which the result 
wid 
nh Z 2 
ee er" le ( . Piece Je 
= 1 
3p, O : 
can easily be obtained. This can be written as 
3/2 
ee Vv 
tS sem 8 
ee 
ca 


A review of the steps taken to arrive at this result shows that 
even had some other dependence on pressure been used, the functional 


dependence does not change because of a temperature effect. 


il 3h 


A theoretical treatment has given several predictions 
as to the operation of a mobility-limited thermionic diode as a 
microphone. The results may be summarized as follows. 
1. The diode d.c. current will vary with pressure and plate 
voltage as 
Wee Se k Sam 2 ae Ss a. 
1 iz Z 
Zee FOXY a Simple circuit lake Figure zi, the Sensitivity es. 


1s given by 


Mg meal DC 
Se V e 
P pe 
1+ k 
2 V-V 
DC 
B. The sensitivity ¢5 s@meeicr coheme wei le fall off “iaeb ape ve- 
length as 
Samet cic, x 
(27h/2 ) y 


2. Experiment 

An experiment was performed using several diodes which 
investigated the first two predictions. The frequency response 
was not tested. 

a. Experimental Arrangement and Procedure 

The diodes used were of similar construction. They 

consisted of a tantalum cylinder, 1 cm long by 1 cm diameter, as 
plate and a tungsten or thoriated tungsten filament as cathode. 
The diodes were mounted several decimeters from the tube ends and 


the sound source, ina tube 4 cm in diameter. A current-modulated 
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glow discharge acted as sound source. The a.c. signal was 
measured using a circuit like Figure 6 with a differential amp- 
lifier to amplify the signal which was then observed with an 
oscilloscope or with a boxcar integrator. 

After outgassing the system and activating the fila- 
ments, space charge limited operation in a high vacuum was verified. 
This was done by measuring the plate current as a function of plate 
voltage to a value of 30 volts, which was twice the highest ex- 
pected operating voltage. Ifa log I, log V plot gave a straight 
line, the selected operating point was considered valid. If not 
and the current dx opped below a straight line at 30 volts, then 


the diode was considered emission limited and appropriate steps 


tem. iiLmthe operating pOine waS Valad, then that filament 
voltage was maintained throughout the calibration. For the re- 


mainder of the calibration, the pressure was changed and the I 
versus V measurement made. For this portion of the calibration, 
1f the log log plot showed that the current at higher voltages was 
above the linear fit to the lower values, than ionization was 
taking place and subsequent operation below these values of plate 
voltage was required. 

After calibration, the a.c. response to various pres- 
sure waveforms was noted. The sensitivity was measured with a 
boxcar integrator set to read the same point in the signal as the 
d.c. Signal was varied by changing the resistance of the load 
resistor. The sound source amplitude was kept constant during the 


measurements. 
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Dee Noise 

The diode was found to be susceptible to electrical 
noise from external and internal sources. The external noise was 
reduced to an insignificant level by shielding with heavy copper 
foil. The foil was wrapped around the tube at its surface, ex- 
tending beyond the diode approximately 5 cm at both ends. The 
Shield was grounded with a separate head. All leads to the diode 
were also shielded. 

In addition to the external noise, the diode was sub- 
ject to an internal noise which seemed due to spurious ionization. 
This noise was characterized by fast rising positive spikes, re- 
presenting a sudden increase in current, followed by a more gradual 
fall. These spikes always appeared when contamination was suspected. 

Ge. D.Csé@maracteristics 
Dayton et.al. [53] reported results for a diode oper- 


ating in N, whereas Argon and Neon were used in this experiment. 


a 
The results obtained were similar to those of Dayton. A typical 
plot for Argon and Neon is shown in Figure 22. The values of Ky 
and Ky obtained are displayed in Table IV along with information on 


Ny taken from Dayton's paper. The values shown are the averages 


computed from data for several diodes for the pressure range 5 to 


PUMmLorr. 
TABLE IV 
Gas Ky Ky k 
Ne .60 + .01 1.67 + .04 1.5 
AX 63 + .02 1.03 + .09 
N 62 1.98 7 


ee 








d. Sensitivity 
A plot of signal voltage versus d.c. voltage for con- 
stant acoustic pressure constitutes a measurement of the sensitivity 
in arbitrary units. Figure 23 shows the results of such a plot. 
The solid curve was calculated using equation (A.17) with kK, and 
the amplitude determined by a least squares fit to the data points. 
The upper half of the points were neglected in the fit because of 
a definite shift in the slope of a linear plot of the data when the 
plate voltage became comparable to the filament voltage. 
3. Discussion of Results 

The Sania of Ky WS Very Similan seneal wen ee gases but 

there is a large difference in the values for ke - Since the 


theeretical tre 


{> 


tment assumes that the pressure dependence arises 
from the similarity parameter E/p, if Ky is the same, then Ky 
should be the same. As this is not the case, and the analysis 
has also shown that the temperature gradient should not affect 
this, then there must be a pressure dependence other than E/p. 
Von Engle C40] points out that ol eeat es Gases. such as 


N have large in-elastic collision cross sections which causes 


Vad 
Hireir mobility to be constant and their drift velocity to vary 
ieinearly with E/p. This is related to the manner in which#€, the 
average fraction of electron energy lost in a collision, increases 
with E/p. Until inelastic collisions start to take place, # re- 


mains constant, afterwards it increases with E/p. The relation 


between the dependence of K on E/p and the dependence of the 
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Figure 23. . Diode Sensitivity in Dependence of d.c. 
Signal. et qonwe<e pressumeuom 7.35 Torr. 
Vertical scale is arbitrary. The ordinate 
is the d.c. component of the signal voltage. 








mobility-limited diode current on plate potential is a complicated 
one and would be difficult to establish quantitatively. However 
there is a correlation between the experimentally determined values 
of kK, and the slope of the log - log E/p curves as figure 63 of 


Von Engle [40]. For Simplicity, k. 1s defined by 


3 
cee 


which is used merely for comparison. The values of kK, determined 
from Von Engle are shown in Table [yfor comparison. 
4. Conclusion 

The mobility-limited diode nicrophone should serve to fill 
a need in plasma-acoustic studies. It is simple and inexpensive 
meeconstruct and requires no complicated circuitry. By virture of 
changing the d.c. signal, sensitivities of different diodes can be 
matched for direct comparison of signals. It is sensitive. 
Sensitivities higher than .27 mV/microbar (360 mV/torr) at 5 torr 
ambient pressure, can easily be attained. It is transparent to 
sound and offers little interference to light in the axial Girection- 
Calibration at zero frequency is simple and easy to perform. 

It is susceptible to noise from contamination. This is a 
drawback which conventional microphones do not have. Its use 
would be limited to inert gases or those where chemical activity 
with the hot tungsten filament would not be expected. Even with 
these drawbacks, a diode microphone can be used for many experi- 


ments where other techniques can not. 
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APPENDIX B 


THE USE OF BUTT SEARS POR Pee Pik @vuGis 

In order to seal electrode supports and feedthroughs into 
glass with a minimum disturbance of the glass surface, a special 
procedure was required. Kovar wire, 7O mil. diameter, was used 
for all supports and feedthroughs. Tantalum sheet 5 mil., was 
used ae the electrodes. Diode filament leads were wrapped with 
fine rhenium wire for ease in spot welding. The electrode supports 
went along the electrode in the axial direction and were spot 
welded in several places. In addition, a small cut in the elec- 
trode was made at the edge and the support wire was bent back and 


crimped to hold the electrode in a mechanical grip. The small 


) 


cut was made to prevent the wire from extending beyond the edge of 


a 


the electrode. The materials were cleaned prior to assembly and 
the components were cleaned again just before sealing in accordance 
with procedures recommended by Rosebury (56). 

For final assembly and sealing, the components of an element 
were assembled on a special.jig. The jig was designed to hold the 
Support wires mee feedthroughs in a plane which would be normal to 
the axis of the glass tubing. After assembly of the components the 
jig was attached to one piece of the glass tubing, earlier mounted 
in the lathe. The jig was placed so that the feedthroughs were 
flush against the edge as shown in Figure 24. The edge used for 
the seals worked best if it was diamond cut. Fire cut edges were 


susceptible to leaks. The second piece of tubing was matched with 
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the edge of the first as shown in Figure 25. A low pressure 

purge of nitrogen was or to prevent the flame from penetrating 
the tube and contaminating the surfaces. After the purge was on 
and the two pieces and jib were turning in the lathe, heat was 
applied and the two ends pushed toward each other. When suffi- 
ciently hot, the glass would flow around the leads and the two 

ends would butt together, forming a butt seal. A small ring of 
7052 was then beaded around the base of the leads, for insurance. 
The complete tube was annealed in an oven following the recommended 
temperature cycle. For the tube used in the experiment, due to the 
nature of the main discharge filaments, annealing was done with 

the tube being evacuated by a diffusion pump. During pre- 
liminary work, several seals were not annealed in the oven, but 


merely flame annealed on the lathe with no later cracking. 
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Figure 24. Diode in Process of Assembly, 
Showing Placement in Jig. 





Figure 25. Diode in Process of Assembly, 
Showing Assembly in Lathe. 
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APPENDIX C 


TAPE RECORDER SPECIFICATIONS 


FM Recor d/Repr oduce 


Input Sensitivity: (High Level) 
(Low Level) 
Input Impedance: (High Level) 


(Low Level) 


1 to 20 volts peak-to-peak 
fo aa 40% deviation 


OZ030 tom volt peak-to-peak 
for + 40% deviation 


20k ohms or greater, unbalanced 
to ground 


50k ohms or greater, unbalanced 
to ground 


Carrier Frequency, Response, Signal-to-Noise Ratio: 


High Level Center Frequency 


73/7) ips 5OCmn2 
B/D) ips 5 khz 
Ege ips 50 khz 


Low Level 


=) ips 500 hz 

e275 ips Sch 

bye9 ips 50 khz 
Gain Drift: 


Center Frequency Drift: 


Linearity: 


Output Level: 


Output Impedance: 


Harmonic Distortion: 
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Response S/N 
DC to 100 hz +1 db 35 db 
DG tod, khz 2 aegdb 38 db 


Deutor 10 khziee dibueeymed 


Sey MOON rans | Clls 315) ele 
Sez tow is kizen, 1edbmeeg omc 
Seniz to TO Maha el clomeaaal 


2% or less of full (temp. range 
Seale Bo oe 1k Oats 
Die. Vase Senile ee ee 
range 110 
a vac to 125 
vac) 


1% of best straight line through 
Zero 


2.8 v p-p for + 40% deviation 
into 1,000 ohms 


50 ohms or less 
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